pHD-fll25  531  THE  USE  OF  TIME  SERIES  ANALVSIS  TO  DEVELOP  SPARES  lA- 

REQUIREMENTS  FORECflSTS(U)  AIR  FORCE  INST  OF  TECH  ^ 

WRIGHT-PRTTERSON  flFB  OH  SCHOOL  OF  SVST.  .  L  D  TflVLOR 
UNCLASSIFIED  SEP  83  AFIT-LSSR-94-83  F/G  12/1  NL 


Wright-Patterson  Air  Force  Base,  Ohio 


THE  USE  OF  TIME  SEMES  ANALYSIS  TO 
DEVELOP  SPARES  REQUIREMENTS  FORECASTS 

Larry  D.  Taylor,  GS-ll 

LSSR  94-83 


'*  ••  i  » 

.J  * 


**  V  / 


KK! VSCSTTTVrrrWT 


ST  V, 


The  contents  of  the  document  are  technically  accurate,  and 
no  sensitive  items,  detrimental  ideas,  or  deleterious 
information  are  contained  therein.  Furthermore,  the  views 
expressed  in  the  document  are  those  of  the  author (s)  and  do 
not  necessarily  reflect  the  views  of  the  School  of  Systems 
and  Logistics,  the  Air  University,  the  Air  Training  Command, 
the  United  States  Air  Force,  or  the  Department  of  Defense. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Doit  Entered)  _ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


t.  REPORT  NUMBER  )2.  GOVT  ACCESSION  NO.j  i-  RECIPIENT’S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


t.  REPORT  NUMBER 

LSSR  94-83 


«.  TITLE  fan d  Submit) 

THE  USE  OF  TIME  SERIES  ANALYSIS  TO 
DEVELOP  SPARES  REQUIREMENTS  FORECASTS 


5.  TYPE  OF  REPORT  b  PERIOO  COVERED 

Masters  Thesis 


«.  performing  org.  report  number 


1.  AUTMORfa) 


Larry  D.  Taylor,  GS-11 


S.  PERFORMING  ORGANIZATION  name  and  AOORESS 


S.  CONTRACT  OR  GRANT  NUMSERfa) 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  b  WORK  UNIT  NUMBERS 


School  of  Systems  and  Logistics 

Air  Force  Institute  of  Technology,  WPAFB  OH 


1 1.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Department  of  Communication 
AFIT/LSH,  WPAFB  OH  45433 


12.  REPORT  DATE 

September  1983 


IS.  NUMBER  OF  PAGES 


4.  MONITORING  AGENCY  NAME  B  AOORESS (II  dllltrtnt  from  Controlling  Oltlet)  I  IS.  SECUPITY  CLASS,  (of  tbit  report) 


ISa.  OECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (a!  tbit  Rtport) 


Approved  for  public  release;  distribution  unlimited 


{  *7.  OtSTftiBUTIOM  STATEMENT  (of  thm  mbttrmct  #n  ffd  In  8toek  20,  It  different  from  Report) 


IB.  SUPPLEMENTARY  notes  «,  ppia  MW  AF3lsiM7 

IXjntE  Wff IAVEH 

Oa«B  tar  Haaaardi  and  Proi^oM,  Ctttitpm,-, 
AU  7 tret  larttute  ol  Tachnoloffr  (ATC) 
.WN<lu-Pattarao«  ATS  03  4M33 


IS.  KEY  WORDS  (Continue  on  rtvtrtt  tldt  II  ntcttttry  and  Identity  by  block  number) 

Demand 

Demand  Rates 
Forecasting 

Requirements  Forecasting 
Time  Series  Analysis 


20.  ABSTRACT  (Continue  on  rtvtrtt  tldt  II  ntcttttry  end  Identity  by  block  number) 

Thesis  Chairman:  Herbert  A.  Stewart,  Major,  USAF 


i  i 


DO  l  jan *71  1473  COITION  OF  I  NOV  SS  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  'When  Dft 


m 


1KCUWITY  CLASSIFICATION  OF  THIS  PAGCfWiwi  Oaf  Bnltfd) 


rhere  is  a  growing  suspicion  that  a  linear  relationship  for  demands 
per  flying  hour,  as  currently  used,  is  not  an  accurate  assumption. 
Fhe  spare  parts  data  used  for  stock  level  maintenance  may  not, 
necessarily,  conform  to  the  general  assumption  of  a  linear  relation¬ 
ship  for  demands  on  those  parts  per  flying  hour.  There  is  a  concern 
that  support  based  on  this  peacetime  assumption  may  fall  short  in 
a  wartime  situation.  Demands  per  flying  hours  are  inherently  a 
factor  of  time,  but  this  factor  is  not  presently  considered  in  comp¬ 
uting  forecasts.  Time  series  methods  are  compared  to  a  model  of 
the  same  data,  produced  without  the  benefit  of  time  series  methods, 
to  demonstrate  that  using  the  f linearity  assumption1!  in  forecasting 
can  lead  to  errors.  L  ’ 


UNCLASSIFIED 


LSSR  94-83 


THE  USE  OF  TIME  SERIES  ANALYSIS 
TO  DEVELOP  SPARES  REQUIREMENTS  FORECASTS 


A  Thisit 

Presented  to  the  Faculty  of  tha  School  of  Systems  and  Logistics 
of  tho  Air  Fores  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the  Requirement  for  the 
Degree  of  Master  of  Science  in  Logistics  Management 


Larry  D.  Taylor,  BS 

es-ii 


September  1983 


Approved  for  public  release! 
distribution  unlimited 


This  thesis,  written  by 


Mr.  Larry  D.  Taylor 

has  been  accepted  by  the  undersigned  on  behalf  of  the 
faculty  of  the  School  of  Systems  and  Logistics  in  partial 
fulfillment  of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  LOGISTICS  MANAGEMENT 

Date:  28  September  1983 


READER 


TABLE  OF  CONTENTS 


Page 


COMMITTEE  APPROVAL  SHEET  .  ii 

LIST  OF  FI6URES .  v 

LIST  OF  TABLES .  vi 

CHAPTER 

1.  INTRODUCTION 


Background . 1 

Resource  Determination . 3 

Problem  Statement . 3 

Research  Objectives . 6 

Justification  .  7 

Research  Question  .  8 

Scope  and  Limitations .  8 

Methodol ogy  .  9 

Preview . 10 

2.  LITERATURE  REVIEW 

THE  SPARES  DETERMINATION  PROCESS  .  11 

Interim  Contractor  Support  .  11 

Initial  Provisioning  . .  12 

Initial  Spares  Support  List  .  13 

Demand  Levels  . .  14 

CURRENT  RESEARCH  .  IS 

Sortie  Length  and  Linearity  .  17 

Inconsistencies  in  Readiness  Goals  .......  20 


Demand  Rates  and  the  Linearity 
Assumption  . . 

TIME  SERIES  ANALYSIS  . 

USAF  Study  . 

Use  of  Time  Series  Analysis  ... 

3.  METHODOLOGY . . 

Time  Series  Analysis  -  Overview  .. 

Prewhitening  the  Univariate  Series 

The  Final  Prewhitened  Univariate 
Model  . 

The  Multivariate  Model  . 

The  Transfer  Function  Model  . 

An  Example  of  a  Transfer  Function 
Using  Raw  Data  . 

4.  CONCLUSIONS  AND  RECOMMENDATIONS  . 

Findings  . 

Problems  with  Using  Time  Series 
Analysis  . 

Concl usi ons  . . 

Recommendat ions  . 

APPENDICIES  . 

A.  THE  DATA  . 

B.  THE  UNIVARIATE  MODELS  . 

C.  THE  MULTIVARIATE  MODELS  . 

SELECTED  BIBLIOGRAPHY  . 

A.  REFERENCES  CITED  . 


B.  RELATED  SOURCES 


LIST  OF  FIGURES 


m 


n 


Figure 


'  1 


Page 


Deviations  from  the  Mean  -  Non-St at ionary  .. 


Deviations  -from  the  Mean  -  Stationary 


Autocorrelation  Function 


Partial  Autocorrelation  Function 


Sumeary  of  Model 


Power  Spectrum 


Best  Fit  Periodograe  of  Actual  Flying  Hours 

Autocorrelation  Function  of  ARIMA  (2,1,2) 
Model  . 

Partial  Autocorrelation  Function  of  ARIMA 
(2,1,2)  Model  . . . 


Summary  of  Actual  Flying  Hours  Model  - 
ARIMA  (2,1,2)  . 


Cross  Correlation  Plot 


Estimated  Impulse  Weights 
Transfer  Function  Summary 


Transfer  Function  Periodogram 


Histogram  for  WUC:14FB0  Transfer  Function  . 


Cross  Correlation  for  WUC: 14FB0 


Deviations  from  the  Mean 


Deviations  From  the  Mean 


Periodogram  for  the  ’Old  Method’ 


Hi stogram 


Summary  of  Model 


Cresc  r'rrr  scions 


CHAPTER  1 


INTRODUCTION 


Background 

An  ever  growing  concern  o-f  the  U.S.  Air  Force  is  the 
question  o-f  whether  our  -forces  are  well  enough  equipped  and 
trained  to  Maintain  the  deterrent  stature  that  the  United 
States  has  enjoyed  -for  so  long.  I-f  the  United  States  should 
have  to  go  to  war,  then  how  well  will  the  peacetime  plan¬ 
ning  support  and  sustain  our  -forces?  Peacetime  planning  and 
training  is  based  on  a  "best  estimate”  o-f  what  we  might  ex¬ 
pect  to  -face  in  the  "next  war.”  The  increasingly  complex 
weapon  systems,  designed  to  counter  the  threat  o-f  potential 
enemies,  make  the  problem  o-f  supporting  and  sustaining  for¬ 
ces  more  acute  than  ever  before  in  our  history.  Faced  with 
the  additional  constraint  of  limited  funds  to  support  a 
weapon  system,  the  Air  Force  must  be  able  to  procure  the 
materials  necessary  to  support  the  forces  in  the  most  effec¬ 
tive  manner  possible. 

The  cost  of  a  weapon  system  is  so  high  that  it  is  vir¬ 
tually  Impossible,  as  well  as  impractical,  to  have  a  “one  of 
everything"  type  of  supply  system  to  support  it.  The  limited 
funds  available  to  the  Department  of  Defense  (DoD)  for  sup¬ 
port  and  maintenance  of  its  weapon  systems,  requires  that 
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prudent  judgement  be  exercised  when  considering  what  and  how 
much  to  buy.  A  large  percentage  of  the  dollars  spent  during 
the  life  of  a  weapon  system  are  for  its  maintenance  and  re¬ 
placement  or  repair  of  parts. 

The  nature  and  purpose  of  a  weapon  system  demand  a  de¬ 
sign  such  that  individual  components  can  be  quickly  removed 
and  replaced  to  minimize  its  down  time.  These  components  are 
often  expensive,  and  constitute  a  large  portion  of  the  cost 
of  the  system.  The  complexity  of  the  components  also  means 
that  maintenance  or  repair  is  costly.  The  Air  Force  has  cat¬ 
egorized  components  into  two  kinds  of  items,  namely,  consum¬ 
ables  and  recoverables. 

Consumable  items  are  either  consumed  in  use  or  otherwise 
lose  their  identity  by  incorporation  into  another  assembly. 
They  are  generally  low  cost  and  most  often  supplied  in  eco¬ 
nomic  order  quantities  (EOQ) .  Recoverable  items  are  general¬ 
ly  high  cost,  retain  their  individual  identity  (usually  via 
serial  numbers),  and  are  repairable  or  reuseable.  This  group 
of  items  is  further  broken  down  into  Line— Replaceable-Units 
(LRU'S),  which  can  be  removed  and  replaced,  as  a  unit,  to 
correct  a  deficiency  or  malfunction  of  a  weapon  system  (or 
other  major  items  such  as  equipment).  These  are  broken  down 
into  Shop-Replaceable-Units  (SRU’s),  or  sub-components  mak¬ 
ing  up  LRU’s.  The  terms  LRU  and  component  will  be  used  syn¬ 
onymously  in  this  thesis. 

The  LRU's  and  SRU's  are  repairable  at  various  levels  of 
maintenance,  i.e. ,  operational  (base),  intermediate  and  de- 


pot.  Base  level  maintenance  is  that  level  of  maintenance 
determined  to  be  within  the  capabilities  of  an  operating 
base  or  wing  with  its  assigned  skills  and  equipments.  Inter¬ 
mediate  level  maintenance  is  normally  accomplished  at  the 
major  command  level.  It  may  be  performed  at  a  fixed  or  mo¬ 
bile  location,  within  the  geographic  area  of  the  command  and 
may  be  collocated  with  an  operating  base.  Depot  level  main¬ 
tenance  is  a  more  specialized  maintenance  which  consoli¬ 
dates  highly  technical  skills  and  specialized  equipments  for 
major  repair  and  modification  of  systems  and  recoverable 
items.  After  the  initial  procurement  of  a  weapon  system,  the 
depots  are  responsible  for  maintaining  the  overall  supply 
and  distribution  of  spare  parts  and  materials  for  logistic 
support  of  assigned  weapon  systems,  or  equipments,  for  their 
useable  life.  This  includes  the  additional  procurement  of 
LRLI’s  and  other  parts  to  maintain  a  determined  level  of 
usage. 

Resource  Determination 

Demands  per  flying  hour  data  is  currently  being  used  by 
the  Air  Force  to  forecast  spares  requirements.  The  develop¬ 
ment  of  complex,  dynamic,  computer  models  to  assist  in  the 
management  of  resources  has  given  rise  to  the  question  of 
how  accurate  forecasts  based  on  this  data  really  are.  Demand 
per  flying  hour  (D/FH)  is  a  statistic  based  on  the  require¬ 
ment  to  replace  some  part  on  an  aircraft  after  a  specified 


period  of  operation.  This  period,  measured  in  flying  hours. 


V»  .  *  .  *  .  •  ,v 


and  th«  demand  ratio  a re  projected  according  to  the  histori¬ 
cal  number  of  replacements  for  each  individual  component. 
Most  computations  of  stock  levels  assume  a  steady  state  with 
constant  demand  and  service  rates  exists.  The  assumption  has 
proven  both  reasonable  and  convenient  during  peacetime 
(11:8). 

When  calculating  the  level  of  stockage  for  an  item,  the 
general  assumption  has  been  a  1:1  linear  relationship  be¬ 
tween  its  demand  and  flying  hours.  That  is,  as  the  histori¬ 
cal  records  of  components  are  compiled,  and  the  data  is  av¬ 
eraged  over  the  cumulative  flying  hours  for  a  fleet  of  air¬ 
craft  using  that  component,  the  resultant  computation  shows 
that  "X"  number  of  components  have  been  requisitioned  for 
every  "X"  number  of  flying  hours.  The  relationship  of  de¬ 
mands  on  an  LRU  as  a  function  of  flying  hours  has  important 
bearing  on  the  stock  level  of  that  LRU. 

The  forecasts  use  mathematical  models  calibrated  with 
data  collected  in  a  peacetime  operational  environment.  For 
the  most  part,  this  D/FH  ratio  has  been  adequate.  This  ratio 
could,  conceivably,  be  less  linear  than  the  relationship  gen¬ 
erally  assumed  for  the  peacetime  weapon  system.  Some  prelim¬ 
inary  investigations  argue  that  breaks  per  flying  hour  will 
be  sub 1 inear  ("aircraft  get  better  the  more  they  are  flown") 
for  some  components  and  super linear  (parts  fail  more  the 
more  they  are  flown)  for  others  (8:6).  These  investigations 


have  examined  'war  games’  or  surge  exercises  (intended  to  be 
as  realistic  as  possible),  for  their  overall  affect  on  sup- 


ply  and  -found  that,  "Overall,  total  demands  per  sortie  were 
comparable  to  (lower  than,  etc.)  normal  peacetime  -flying 


(8:6)."  A  surge  exercise  is  designed  to  demonstrate  a  unit’s 
ability  to  sustain  the  high  level  o-f  flying  activity  that 
may  be  expected  in  war.  It  seems  intuitive  that  a  wartime 
environment  should  place  different  levels  of  demands  on  the 
inventory  system  than  peacetime. 

Problem  Statement 

Many  maintenance  personnel  who  have  worked  with  spe¬ 
cific  systems  for  some  time  believe  that  as  their  aircraft 
get  more  and  steady  usage  that  the  rate  of  failure  for  com¬ 
ponents  is  less,  or  sometimes  more,  frequent  than  the  supply 
system  provisions  for  (8:6).  If  in  fact  the  D/FH  ratio  is 
other  than  1:1,  a  more  accurate  forecasting  method  should 
provide  dollar  savings,  as  well  as  enhance  the  availability 
and  readiness  of  aircraft. 

The  spare  parts  data  used  for  stock  level  maintenance  is 
based  on  a  steady  state,  constant  demand,  peacetime  environ¬ 
ment,  which  lends  itself  to  the  present  usage  of  a  1:1  line¬ 
ar  ratio  for  the  D/FH  for  stock  levels.  The  data  is  made 
linear  by  using  the  equation  Y*  -  a^+a^+e*  which  is  essen¬ 
tially  the  same  as  that  for  a  straight  line.  Here  Y«  is  the 
demand  at  period  t,  the  ai’s  are  the  values  of  the  observa¬ 
tions  at  X«,  which  is  the  value  of  the  flying  hours  at  peri¬ 
od  t,  e«  is  the  error  factor  (17).  (This  type  of  notation  is 
used  throughout  the  thesis).  The  demand  data  is  accumulated 


and  is  than  dividad  into  tha  cumulative  flying  hours  which 


gives  a  li  near,  1:1  ratio.  With  the  emphasis  on  the  more 
realistic,  dynamic  world,  this  “linearity"  may  not  hold 
true.  The  data  is  time  related,  and  the  calculations  used  in 
the  stock  leveling  do  not  take  this  into  consideration. 

There  are  many  computerized  models  and  forecasting  meth¬ 
ods  available  now  that  theoretically  should  provide  a  more 
accurate  forecast  for  spares.  This  thesis  proposes  to  show 
that  time  series  analysis  is  a  more  accurate  method  than 
using  the  linear  assumption.  Time  series  analysis  is  a  pow¬ 
erful  forecasting  method  that  takes  into  consideration,  be¬ 
sides  the  time  factor,  such  things  as  seasonal  fluctuation 
in  demands.  It  is  able  to  analyze  data  in  the  manner  that  it 
actually  occurs  so  the  forecast  figure  will  account  for  the 
fluctuations,  seasonality  and  other  time  related  facts  that 
Influence  a  forecast  of  future  requirements. 

Bwtrch  QfritcUvft 

This  study  will  use  time  series  analysis  with  selected 
F-16  LRU  data  and  demonstrate  a  method  of  forecasting  future 
requirements  for  those  parts.  The  method  is  applicable  to 
any  similar  time  related  data.  An  additional  outcome  of  the 
study  is  the  intent  to  demonstrate  that  the  number  of  de¬ 
mands  per  flying  hour  is  not  a  linear  ratio.  Therefore,  the 
present  methods  of  forecasting  requirements  may  not  provide 
the  most  cost  effective  forecast,  nor  provide  the  Air  Force 
with  an  optimal  spares  stockage. 


"The  objectives  and  assumptions  employed  in  every  stage 
of  the  resource  allocation  process  ought  to  be  consistently 


derived  from  system-level,  out put -oriented  goals  (9:7). “ 

This  study  supports  the  effort  to  better  manage  spares 
and  repair  parts;  and  therefore,  better  support  the  using 
activities.  The  premise  of  this  thesis  is  that  spare  parts 
required  to  maintain  an  aircraft  may  not,  necessarily,  con¬ 
form  to  the  general  assumption  of  a  linear  relationship  for 
demands  on  those  parts  per  flying  hour.  Spare  parts,  or  the 
lack  thereof,  has  always  been  an  issue  of  concern,  causing 
the  problem  of  allocating  shortages,  rather  than  resources, 
because  requirements  inevitably  exceed  funds  availability. 

An  improved  forecast  might  alleviate  shortages,  improve 
availability  of  ready  aircraft,  and  allow  funds  to  be  spent 
in  a  more  efficient  manner.  The  demands  for  parts  and  the 
flying  hours  are  inherently  a  factor  of  time,  but  this  fac¬ 
tor  is  not  presently  considered  in  computing  forecasts.  Re¬ 
quirements  inevitably  grow  over  time,  while  at  the  same  time 
funds  availability  tends  to  decline,  which  puts  pressure  on 
the  managers  to  delete  unjustified  requirements  (9:9).  All 
of  this  lends  itself  to  the  need  for  more  accurate  forecast¬ 
ing,  so  that  the  required  parts  are  available  while  funds 
are  judiciously  spent. 
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There  is  a  small  amount  of  literature  and  a  growing  sus¬ 


picion  that  a  linear  relationship  (demands  per  flying  hour), 
as  currently  used,  is  not  an  accurate  assumption.  Some  cur¬ 
rent  research  of  the  Rand  Corporation  has  clearly  demonstra¬ 
ted  that,  at  least  for  some  LRU’s,  the  demands  per  flying 
hour  is  not  a  linear  ratio.  Time  series  analysis,  on  the 
other  hand,  does  not  rely  on  cumulative  averages,  but  uses 
the  actual  data  as  it  occurred  through  time  for  making  fore¬ 
casts.  This  thesis  responds  to  the  question  of  whether  a 
better  method  to  forecast  spares  exists,  and  if  so,  how  much 
improvement  can  be  achieved  in  forecasting?  In  attempting  to 
answer  this  question,  the  method  to  arrive  at  a  correct  ra¬ 
tio  for  the  demands  per  flying  hour  to  use  in  forecasting 
future  requirements  will  be  demonstrated.  The  method  can  be 
applied  to  any  aircraft  or  system  that  has  accumulated  data 
over  time. 


Scocit  *nd ,  i  roe. 

To  accomplish  the  study,  five  LRU’s,  which  were  among 
the  "top  ten"  items  impairing  the  F-16  mission  capability 
(MICAP)  for  the  week  of  12  June  1983,  will  be  used  as  the 
data  base.  The  fact  that  these  LRU’s  are  MICAP  items  implies 
that,  for  what  ever  reason,  the  original  spares  forecast 
proved  inaccurate.  These  five  LRU’s  provide  ample  data  to 
demonstrate  the  method  of  time  series  forecasting  spares  re— 


quirements.  Appendix  A  contains  a  listing  of  the  LRU's  by 


work  unit  cods  (WUC) ,  actual  -Flying  hours  (by  month,  06/79 
to  04/83  inclusive),  and  the  removals  or  demands  -For  the 
LRU's  each  month.  This  is  the  same  type  o-F  data  that  would 
be  accumulated  and  averaged  to  get  the  linear  ratio  -For  the 
current  -Forecasting  method. 


Hethodol oav 

The  Box— Jenkins  method  o-F  time  series  analysis  will  be 
used  -For  this  thesis.  Time  series  analysis  o-F  data  which 
include  an  independent  variable  (flying  hours)  and  a  depen¬ 
dent  variable  (demand)  is  accomplished  using  a  transfer 
function.  Each  variable  is  first  prewhitened  by  fitting  an 
appropriate  model  to  the  data.  Prewhitening  is  done  to  clar¬ 


ify  the  relationship  between  flying  hours  and  demand. 


For  a  given  pair  of  series,  there  are  interrela¬ 
tionships  between  the  series  and  intrarelationships 
within  the  individual  series.  The  intrarelationships 
are  summarized  in  the  autocorrelation  structure  of  the 
individual  series  CIS: 3783 


The  prewhitened  model,  then,  has  all  the  autocorrelation 
within  the  series  removed  and  only  "white  noise"  remains. 
White  noise  is  analogous  to  a  series  that  is  completely  ran¬ 
dom,  i.e. ,  there  is  no  pattern  whatsoever  contained  in  the 
series.  This  is  necessary  before  the  interrelationship,  or 
cross  correlation,  between  the  variables  can  be  found.  After 
a  prewhitened  model  is  accomplished  for  both  the  variables 
they  are  combined  into  a  single  transfer  function  model  to 


CHAPTER  2 


LITERATURE  REVIEW 


This  rev lew  starts  with  s  look  at  how  spars  parts  ars 
dster mined  for  ths  support  of  a  wsapon  system,  and  will  bs 
followsd  by  a  rsview  of  ths  current  litsraturs  docunsnting 
rsssarch  on  sons  of  ths  pr obi sms  sncountsrsd  providing  this 
support.  Finally  a  ssction  on  tins  ssriss  analysis  and  why 
it  should  bs  considsrsd  as  part  of  ths  solution  to  ths  sup¬ 
port  problems  will  bs  presented. 

THE  SPARES  DETERMINATION  PROCESS 

Contractor  Support 

During  ths  development  phase  of  a  wsapon  system,  design 
modifications  and  changes  can  occur  frequently  until  the 
updated  requirements  for  the  system  are  met.  As  a  part  of 
the  contract,  the  contractor  is  usually  assigned  Interim 
Contractor  Support  <ICS)  during  the  demonstration  and 
validation  (D&V)  phase  and  until  the  Air  Force  is  able  to 
organically  support  the  weapon  system  (2:1).  The  Air  Force 
Program  Manager  (PM)  and  the  Deputy  Program  Manager  for 
Logistics  (DPML)  analyze  the  risk  of  the  components  (i.e. , 
to  defer  total  investment  until  after  the  D&V  phase  or  until 


the  design  is  stabilized),  to  determine  which  components 
should  be  assigned  ICS.  This  list  of  ICS  components  is 
further  refined  during  the  full-scale  engineering 
development  and  contracted  for  during  the  production  and 
development  phase  (2:9). 

Initial  Provisioning 

The  initial  provisioning  of  parts  for  a  new  weapon  sys¬ 
tem  is  developed  after  the  contract  to  build  the  system  has 
been  awarded.  Initial  provisioning  is  defined  in  AFLCR  65-5 
as  a  "management  process  for  determining  and  acquiring  the 
range  and  quantity  of  support  items  necessary  to  operate  and 
maintain  and  end  item  of  material  for  an  initial  period  of 
service  (1:1-1)."  The  initial  period  is  defined  "to  begin 
with  the  preliminary  operational  capability  of  the  weapon 
system  through  the  lead  time  required  for  the  component, 
plus  three  months,  or  a  12  month  minimum  (1:1-1)."  This  ini¬ 
tial  provisioning  is  designed  to  give  adequate  parts  support 
until  normal  resupply  is  established.  During  this  time  the 
contractor  and  the  Air  Force  jointly  determine  the  items 
that  have  to  be  provisioned. 

All  parts  of  the  weapon  system  are  reviewed  for  current 
availability  in  the  Air  Force  inventory.  If  items  are  avail¬ 
able,  the  Item  Manager  (1M)  is  informed  so  that  stock  levels 
can  be  adjusted  upward  for  the  additional  requirements. 
Items  not  currently  in  Air  Force  inventory  are  assigned  Na¬ 


tional  Stock  Numbers  (NSN)  and  assigned  to  an  IM  for  pro- 


cur ■man t  and  subsequent  management.  Where  an  item  is  already 
in  the  inventory,  the  present  demand  data  is  already  estab¬ 
lished  and  only  needs  adjustment  based  on  initial  -forecasts. 

The  items  that  are  new  to  the  Air  Force  do  not  have  an 
established  demand  and  must  be  initially  provisioned  so  they 
are  available  -for  use  Mhen  the  weapon  system  becomes  opera¬ 
tional.  Consumable  (EOQ)  items  are  provisioned  based  on  es¬ 
timates  using  demand  information  from  similar  items  current¬ 
ly  in  the  inventory.  The  recoverable  items  not  currently  in 
the  inventory  are  provisioned  based  on  a  combination  of  fac¬ 
tors  gained  from  the  contractor  and  factors  developed  by  an 
Air  Force  equipment  specialist.  The  contractor  furnishes 
such  data  as  the  items  mean  time  between  failures  (MTBF) , 
the  mean  time  to  repair  (MTTR),  the  replenishment  cycle  time 
(RCT)  and,  the  lead  time  for  procurement  (3).  The  equipment 
specialist  uses  this  information  to  develop  the  initial 
maintenance  requirements.  The  demand  rate  for  the  item  is 
initially  established  through  the  compilation  of  this  data. 

Initial  Soares  Support  List 

The  Initial  Spares  Support  List  (ISSL)  is  a  listing  of 
recoverable  and  consumable  line  items  used  in  the  repair  of 
the  recoverables.  It  is  intended  to  include  the  parts  neces¬ 
sary  for  base  level  maintenance.  The  ISSL  is  developed  by 
the  depot,  which  is  an  Air  Logistics  Center  (ALC) .  It  is 


developed  during  the  provisioning  phase  and  is  based  on  in¬ 
puts  from  the  contractor,  the  IN,  the  System  Manager  (SM) 


.  .*  -  V  -•  ^ 


and,  the  using  command  that  Mill  "own"  the  weapon  system 
(14).  The  quantity  o-f  recoverable  items  on  the  ISSL  is  the 
additive  -factor  o*f  the  base  order  and  shipping  time  and  the 
base  repair  cycle  time  quantity  (19).  Consumable  quantities 
are  based  only  on  the  base  order  and  shipping  time  quanti¬ 
ties.  The  ISSL  quantities  are  put  into  the  Standard  Base 
Supply  System  (SBSS)  computer  as  the  initial  demand  level. 
This  demand  level,  a-fter  two  years,  is  adjusted  based  on 
actual  demands  -for  the  parts. 

Demand  Levels 

The  demand  levels  -for  recoverable  items  are  determined 
by  the  SBSS,  by  adding  the  base  repair  cycle  quantity,  the 
order  and  shipping  time  quantity,  the  quantity  not  repaira¬ 
ble  this  station  (NRTS)  condemnation  rate,  plus  a  safety 
stock  level.  To  this  figure  a  +0.5  adjustment  factor  is  add¬ 
ed  for  items  costing  S750. 00  or  more,  and  a  +0.9  adjustment 
factor  for  those  costing  less  (1).  The  base  level  demands 
are  forwarded  to  the  ALC  responsible  for  the  weapon  system. 
There,  each  independent  asset  requirement  is  consolidated 
into  the  Recoverable  Consumption  Item  Requirements  System, 
the  D041.  Inputs  to  the  D041  include  major  command  (MAJCOli) 
flying  hour  projections  for  five  years,  and  inputs  from  each 
base  SBSS  which  gives  the  recorded  failures  to  date.  The 
failure  rates  are  applied  to  total  flying  hours  to  get  a 


historical  failure  per  flying  hour  rate.  The  failures  are 
then  applied  to  the  quarterly  flying  hour  projections  to 


obtain  an  estimated  number  of  item  failures  per  future  quar — 
ter.  The  item’s  NRTS  condemnation  rate  is  subtracted  from 
the  projected  number  of  item  failures  to  generate  an  item- 
buy  recommendation.  Any  remaining  failures  are  coded  as  re¬ 
pairs.  The  provisioning  and  War  Readiness  Materials  (WRM) 
requirements  are  input  from  other  automated  systems.  WRM  is 
defined  as  those  assets,  in  addition  to  normal  peacetime 
operating  stocks  (POS) ,  needed  at  the  start  of  a  conflict  to 
support  active  forces  until  industrial  production  can  sus¬ 
tain  combat  requirements  (4:1).  All  other  sources  of  re¬ 
quirements  are  manually  input  by  the  IM.  These  include  For¬ 
eign  Military  Sales  (FMS) ,  retrofit  requirements,  non  pro¬ 
grammed  requirements,  training  and  testing  requirements  and 
the  1SSL. 

CURRENT  RESEARCH 

This  issue  of  "linearity"  led  to  this  study  to  define  a 
better  model  to  use  in  forecasting  future  requirements.  A 
correct  time  series  model  will  eliminate  some  of  the  causes 
of  linearity.  This  thesis  is  mainly  concerned  with  the  meth¬ 
odology  of  using  time  series  methods  as  a  better  tool  in 
forecasting  demands,  or  removals,  for  spare  parts.  Although 
there  does  not  appear  to  be  literature  that  specifically 
addresses  this  issue,  there  is  a  limited  amount  of  litera¬ 
ture  that  does  address  the  linearity  issue.  The  linearity 


issue  and  this  study  are  both  interrelated  in  that  the  as¬ 
sumption  of  linearity,  no  matter  what  the  form,  is  the  cause 
for  many  of  the  logistic  support  problems  being  encountered 
today.  It  is  something  that  has  attracted  the  attention  of 
researchers  who  are  looking  both  for  reasons  why  material 
support  is  not  as  good  as  it  might  be  and  for  better  ways  to 
maintain  the  weapon  systems. 

The  subject  of  ’linearity’  is  a  complex  one.  Only  rela¬ 
tively  recently  in  the  research  is  it  specifically  addres¬ 
sed.  There  are  several  references  that  indirectly  refer  to 
it  and  some  recent  unpublished  Rand  Corporation  studies  are 
more  specific  about  how  the  linearity  of  demands  per  flying 
hour  could  play  a  part  in  spares  computation. 

None  of  the  literature  reviewed  discusses  the  time  fac¬ 
tor,  which  is  an  integral  part  of  the  data  collected  fcr  Air 
Force  requirements,  because  it  is  collected  over  time.  The 
literature  does  discuss  the  interrelationships  of  how  the 
resource  requirements  system  as  a  whole  is  affected  by  the 
nebulous  assumption  of  linearity.  However,  this  thesis  is 
more  concerned  with  finding  an  accurate  model  to  forecast 
time  related  data,  than  with  the  actual  degree  of  "linear — 
ity"  as  it  is  discussed  in  the  literature. 

The  findings  of  the  researchers  demonstrate  that  the 
ratio  of  demands  per  flying  hour  do  change.  They  also  show 
that  the  ratio  is  not  a  linear  one,  for  the  components  in¬ 
vestigated  (7,8,9,10,11,12,16).  Embry  and  Crawford  make  spe¬ 
cific  reference  to  the  issue  of  "linearity"  and  it’s  rela- 


tion  to  the  overall  Air  Force  requirements.  They  define  the 
more  general  type  of  problems  the  issue  of  linearity  poses. 


Peacetime  experience  is  commonly  used  to  develop 
component  demand  rates  (expressed  in  removals  per  fly¬ 
ing  hour)  for  use  in  peacetime  spares  requirements 
computations.  This  same  rate  is  used  as  an  input  to 
Martime  resource  requirements  and  capability  assess¬ 
ment  models  based  on  the  assumption  that  removals  are 
proportional  to  flying  activity.  Since  both  require¬ 
ments  and  capability  estimates  are  sensitive  to  the 
demand  rate  parameter,  they  are  also  sensitive  to  this 
"linearity"  assumption  C10:iii3. 


Sortie  Length  and  Linearity 

Embry  and  Cranford  (10)  examined  a  number  of  past 
studies  that  indicate  the  "linearity"  assumption  is  tenuous 
when  applied  across  operations  having  different  sortie 
lengths.  Also,  some  studies  suggested  the  assumption  may  not 
be  valid  across  periods  with  different  sortie  rates.  Their 
key  point  was  that  if  the  “linearity"  assumption  is  not 
valid,  capability  assessments  that  rely  on  it  must  be 
inaccurate,  and  that  requirements  statements  that  rely  on 
this  assumption  may  produce  the  wrong  mix  of  resources  to 
support  wartime  operations. 

Both  requirements  computation  and  capability  as¬ 
sessment  require  assumptions  to  support  extrapolation 
of  peacetime  experience  to  wartime  levels  and  patterns 
of  activity.  Most  of  the  logistics  models  that  are 
currently  being  used  for  developing  resource  require¬ 
ments  and  assessing  capability  assume  component  de¬ 
mands  are  a  function  of  flying  hours.  Since  both  types 
of  estimates  are  particularly  sensitive  to  this  demand 
assumption,  it  is  important  that  demand  rate  estimates 
be  accurate.  The  "linearity”  assumption,  which  holds 
that  doubling  the  level  of  flying  activity  will  double 
the  demand  for  components,  has  significant  resource 
implications  and  is  a  potential  source  of  bias  in  cap¬ 
ability  projections  C5:ll. 
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The  requirements  computation  -For  spares,  as  described 
above,  considers  the  reported  removals  of  an  item  as  part  of 
the  calculation.  Removals  happen  -For  a  myriad  of  reasons, 
failure  being  the  most  important.  It  is  often  the  case  that 
the  failure  cannot  be  duplicated  (CND)  in  the  shop  environ¬ 
ment,  but  the  removal  is  still  recorded.  Embry  and  Crawford 
point  out  that,  "It  is  important  to  note  that  peacetime  de¬ 
mand  rate  measures  removals,  not  failures  (10:3)."  However, 
the  removal  and  demand  rate  are  considered  synonymous  when 
calculating  spares  requirements.  The  methodology  chapter  of 
this  thesis  will  use  removal  data  for  selected  F-16  LRU's  to 
demonstrate  the  forecasting  technique.  This  will  not  provide 
the  total  solution  to  the  linearity  problem  because  the  as¬ 
sumption  that  removals  are  proportional  to  flying  hours  is 
not  an  accurate  one.  The  time  series  model  will  demonstrate 
this,  and  demonstrate  its  utility  by  showing  a  way  to  "get 
around"  this  linear  assumption. 

Embry  and  Crawford  make  a  point  where  this  assumption 
has  been  used  that,  "If  expected  wartime  removals  are  not 
strictly  a  function  of  the  peacetime  demand  rate  and  wartime 
flying  hours,  both  requirements  estimates  and  capability 
assessments  are  affected  (10:4)".  The  best  solution  then, 
would  appear  to  be  that  "...of  improving  the  basis  for  both 
types  of  estimates  (10:4)" 


Eight  different  studies  that  deal  with  linearity  were 
discussed  in  the  paper  by  Embry  and  Crawford.  In  all  of  the 


studies  the  changes  in  sortie  length  and  rate  were  the  most 
important  variable  in  the  observed  changes  in  demand  rate. 
When  comparing  the  demands  during  the  SALTY  ROOSTER  exer — 
cise  at  Hahn  Air  Force  Base  with  the  rates  during  the  pre¬ 
vious  year,  Crawford  -found  significant  changes  in  a  number 
of  components.  The  removal  rates  ratio  varied  from  .33  to 
2.7,  with  the  higher  rates  belonging  to  items  that  are  not 
exercised  extensively  in  peacetime  but  stressed  more  under 
exercise  conditions.  Other  demand  ratios  dropped  to  a  third 
of  the  peacetime  ratios  (10:8).  In  another  test  of  F-16 
data,  the  ratio  of  surge  to  control  period  malfunctions  per 
sortie  on  the  Fire  Control  and  Weapons  Delivery  System  were 
examined.  Once  again  they  found  that  the  ratio  varied  from 
-.3  to  1.6  when  the  sortie  rate  went  from  .84  to  3.44  sor — 
ties  per  day.  The  inertial  navigation  system  (INS)  experi¬ 
enced  a  1.6  ratio  of  removals  and  all  other  dropped  below 
the  control  average  (10:14).  They  concluded  that  their  find¬ 
ings  reinforced  the  earlier  research,  suggesting  that  chang¬ 
ing  utilization  of  the  aircraft  affects  the  demand  rates 
(10:20). 

All  of  the  research  indicates  that  there  are  problems 
with  using  the  assumption  of  linearity  to  calculate  the  de¬ 
mand  rates  needed  for  forecasting  requirements.  Because  the 
demand  ratio  is  found  to  vary,  rather  than  remain  at  a  fixed 
ratio,  the  assumption  becomes  expensive  in  terms  of  cost  and 
parts  shortages.  It  also  indicates  that  the  mix  of  support 
parts  intended  to  support  wartime  activities  may  not  include 


the  parts  that  will  be  required,  in  the  quantities  required. 
The  spares  kits  may  be  over  stocked  with  some  items  and  un¬ 
der  stocked  with  others}  and  if  the  research  proves  correct 
the  more  expensive,  long  lead  time  items  are  those  that  will 
be  short. 


nconsistencies  in  Readiness  Goals 


In  September  1982,  Embry  (8)  presented  a  briefing  to  the 
DoD  Material  Readiness  and  Sustainability  Symposium  in  which 
he  stressed,  "...  the  need  for  analyses  and  management  ac¬ 
tions  that  recognize  the  interdependencies  of  various  ele¬ 
ments  of  the  aviation  logistics  support  system  (9:iii).M  He 
stressed  that  the  outputs  of  these  interdependencies  be  tied 
to  aircraft  oriented  wartime  performance  measures.  The  man¬ 
agement  and  analyses  actions  that  concern  the  resource  read¬ 
iness  issue  involve  program  adequacy,  the  effects  of  alter¬ 
nate  funding  levels,  and  feedback  and  control.  The  feedback 
and  control  phase  is  where  Embry's  concern  was  stressed 
most.  "The  problem  during  this  phase  is  generally  more  one 
of  allocating  shortages,  not  resources,  because  requirements 
Inevitably  exceed  funds  availability,  particularly  when  re¬ 
quirements  change  (9:6). "  This  problem  of  requirements 
change  is  due,  in  part,  to  fluctuations  in  the  demand  rates 
which  were  discussed  earlier.  Changes  in  requirements  after 
the  determination  of  the  funding  level  make  attainment  of 
readiness  goals  nearly  Impossible. 


These  goals  are  used  to  guide  management  as  well  as  size 
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r •source  requirements.  The  failure  to  meet  the  goals,  is 
where  the  challenge  to  logistic  managers  is  the  greatest. 
Embry’s  comment  about  some  of  the  reasons  for  this  failure 
is  appropriate  here: 

There  are  no  clear  links  between  intermediate  per — 
formance  measures,  such  as  supply  fill  rates  or  main¬ 
tenance  manpower  utilization,  and  wartime  sortie  gen¬ 
eration  capability.  Yet  these  functional  goals  are  the 
basis  for  resource  requirements  computations.  Because 
they  are  not  linked  to  output  measures,  inconsisten¬ 
cies  in  the  assumptions  used  to  size  various  function¬ 
al  requirements  are  not  only  possible,  but  unfortu¬ 
nately  fairly  common  £9:6-73. 

The  key  point  that  he  makes  is  that,  "...the  objectives 
and  assumptions  employed  in  every  stage  of  the  resource  al¬ 
location  process  ought  to  be  consistently  derived  from  sys¬ 
tem-level,  output -oriented  goals  (9:7)."  He  defines  the  real 
outputs  as  those  that,  "...are  reflected  in  weapon  system 
readiness,  not  supply  demand  satisfaction  or  maintenance 
production  rates  (9:9). "  The  ideal  demand  forecast  would 
allow  supply  demands  to  equal  the  requirements  that  are  re¬ 
flected  in  weapon  system  readiness.  However,  because  the 
resource  requirements  must  be  identified  so  far  in  advance, 
to  meet  budgeting  requirements  it  becomes  difficult  to  fore¬ 
cast  accurately,  the  proper  mix  of  items  to  buy.  This  is 
even  more  apparent  when  the  demand  for  those  items  is  not 
the  same  as  was  assumed  when  the  forecast  was  made.  When  the 
forecast  does  not  meet  the  weapon  system  readiness  require¬ 
ments,  readiness  suffers. 


Kanins  and  Craw-ford  (12)  ds-finsd  ths  relationship  be¬ 


tween  linearity  and  demand  rates,  and  how  they  are  used. 

Peacetime  experience  is  commonly  used  to  develop 
component  demand  rates  (expressed  as  removals  per  -fly¬ 
ing  hour)  -for  use  in  peacetime  spares  requirements 
computations.  This  same  demand  rate  is  used  as  an  in¬ 
put  to  wartime  resource  requirements  and  capability 
assessment  models,  based  on  the  assumption  that  remov¬ 
als  are  proportional  to  -flying  activity.  Because  both 
requirements  and  capability  estimates  are  sensitive  to 
the  demand  rate  parameter,  they  are  also  sensitive  to 
this  "linearity"  assumption  C12:iiil. 

1-f  this  "linearity”  assumption  is  not  valid, 
capability  assessments  that  rely  on  it  are  inaccurate. 
Perhaps  more  important,  requirements  computations  that 
rely  on  this  assumption  may  produce  the  wrong  mix  o-f 
resources  to  support  wartime  operations 
C12:iiil. 


Rand  studies  done  prior  to  the  Kami  ns  and  Craw-ford  study 
have  documented  decreases  in  maintenance  actions  per  -flying 
hour  (MAPFH)  that  seem  to  occur  whenever  aircraft  fly  longer 
sorties.  But,  they  point  out,  these  apparent  decreases  are 
ignored  when  computing  war  readiness  requirements.  One  of 
the  reasons  these  decreases  are  ignored  is  the  assumption 
that  MAPFH  will  double  if  flying  doubles.  They  also  point 
out  that  this  generally  conservative  assumption  results  in 
War  Reserve  Material  (WRM)  kits  that  give  more  capability 
than  would  kits  based  on  contrasting  studies  (12:iv) 

Because  funding  for  spares  is  austere,  the  component 
buys  often  result  in  less  of  the  more  expensive  LRU's  and 
more  of  less  expensive  repair  parts.  So  it  is  necessary  that 
the  correct  requirements  are  identified  initially.  This  is 
noted  by  Kami ns  and  Crawford: 


The  argument  that  the  linearity  assumption  results 
in  a  WRM  kit  with  more  capability  may  be  -Fallacious: 
When  maintenance  -Factors  require  large  numbers  o-F  ex¬ 
pensive  LRU’s,  policy  decisions  mitigate  the  dollar 
requirements  so  that  the  kits  will  be  affordable 
C12:iv3. 

No  matter  how  expensive  a  WRM  kit  is  it  must  have  the 
proper  mix  o-F  parts  or  it  will  not  give  the  desired  results 
in  terms  of  capability  .  "Arbitrary  policy  decisions  design¬ 
ed  to  reduce  high  WRM  kit  costs  are  likely  to  make  the  mix 
even  worse  <12:v>."  Because  o-F  this  fact  the  issue  of  line¬ 
arity  becomes  a  central  focus.  Whether  failures  in  aircraft 
systems  (particularly  those  requiring  parts  replacements) 
remain  constant,  when  the  environment  changes  from  peace  to 
war ,  is  the  issue  that  must  be  dealt  with. 

Kami ns  and  Crawford  cite  a  Boeing  study  which  analyzed 
maintenance  data  on  B-52D  operations  in  Southeast  Asia.  The 
findings  of  that  study  found  that  "...the  landing  gear  sys¬ 
tem  sub-linearity  was  so  pronounced  as  to  decrease  total 
failures  when  sortie  rates  were  increased  (12: 1) . "This  sup¬ 
ports  the  findings  of  studies  mentioned  earlier  that  the 
linearity  assumption  may  lead  to  incorrect  buys  and  support 
will  become  faulty.  Another  study  cited,  by  the  Rome  Air 
Development  Center,  showed  that  higher  levels  of  flying  ac¬ 
tivity  were  associated  with  measurable  increases  in  the  ob¬ 
served  MTBF  for  all  types  of  aircraft  studied,  including 
tactical  (12:2). 

At  the  Third  Logistics  Capability  Assessment  Symposium, 
in  March  1983,  Crawford  (7)  reviewed  the  excessive  variation 


in  empirical  demand  rates  -for  some  important  aircraft  compon¬ 
ents.  This  is  a  "...poorly  understood  phenomenon  that  is 
potentially  very  important  to  logistics  capability  assess¬ 
ment  and  requirements  determination  <7:iii>."  He  -first  de¬ 
scribed  the  Poisson  arrival  process,  and  then  related  that 
to  the  way  requirements  and  capability  assessment  models  use 
them.  These  models  make  the  assumption  that  the  time  between 
-failures  is  an  exponential  random  variable  and  within  this 
binomial  approximation  the  sorties  are  independent. 

Crawford  pointed  out  that  the  observed  variations  in  the 
LRU’s  he  checked  are  not  consistent  with  a  Simple  Poisson 
arrival  process.  The  implications  pointed  out  by  Crawford 
are  as  follows: 

There  are  times  when  removals  seem  to  occur  in 
clusters;  and  because  breaks  are  likely  to  occur  in 
clusters  we  must  conclude  that  maintenance  is  not  a 
"memoryless"  process.  But  a  lack  of  memory  is  implied 
by  the  assumption  in  our  models  that  different  sorties 
or  different  time  intervals  are  independent.  Even  more 
troublesome,  sometimes  we  don't  have  any  idea  what  the 
clock  is,  all  we  know  is  that  something  other  than 
sorties  or  flying  hours  seems  to  be  driving  removals 
C7:203. 

It  seems  clear  that  the  methods  the  Air  Force  uses  to 
forecast  spare  parts  requirements  makes  many  assumptions 
that  are  probably  inaccurate.  The  resources  can  be  better 
allocated,  if  the  methods  used  take  into  consideration  the 
more  complex  factors  that  are  present  in  the  actual  data. 
The  new  computerized  data  systems,  such  as  the  Central  Data 
System  for  the  F-16  and  the  Automated  Maintenance  Systen  for 
the  C-S,  are  providing  better  sources  of  data  for  analysis 


than  has  been  available,  so  it  is  becoming  possible  to  take 


these  -factors  into  consideration.  Craw-ford  makes  the  suc¬ 
cinct  observation  that  i-f  we  want  to  "...intelligently  buy 
spare  parts,  even  -for  peacetime,  we  need  to  know  how  much  o-f 
the  variation  is  manageable  and  then  buy  resources  to  cover 
the  rest  (7:20)."  He  goes  on  to  state  how  much  harder  it 
will  be  to  do  this  -for  wartime  support,  especially  when  the 
rami -fi  cat  ions  o-f  that  kind  o-f  environment  are  considered. 

TIME  SERIES  ANALYSIS 

This  section  reviews  a  time  series  analysis  that  was 
recently  completed  on  USAF  C-141  data.  The  methodologies 
used  in  the  article  reviewed  are  generally  the  same  as  that 
used  -for  the  thesis  research,on  the  F-16  data  in  Appendix  A. 
Because  o-f  the  complicated  description  o-f  time  series  analy¬ 
sis,  detailed  methodologies  used  in  these  studies  will  be 
held  over  to  chapter  3  o-f  this  thesis,  where  examples  and 
graphs  will  aid  in  describing  these  methods,  thus  making  the 
description  more  meaningful  and  avoiding  duplication  of 
effort. 


USflF  atucfy 

In  December  1982,  Singpurwalla  and  Talbott  (16)  used 
time  series  analysis  to  investigate  the  interrelationships 
between  alert  availability  and  flying  hours  for  the  C-141 
aircraft  of  the  U.S.  Air  Force.  Their  data  included  a  major 


change  in  operating  hours  which  was  brought  about  by  the 
* 

introduction  of  the  Reliability  Centered  Maintenance  (RCM) 
concept  in  1976.  Their  data  was  also  “contaminated "  because 
of  anomalies  in  the  reporting  procedures  during  the  accumu¬ 
lation  of  the  data  (16:1).  The  technique  that  they  used  was 
a  combination  of  transfer  function  modeling  and  intervention 
analysis. 

In  a  previous  study,  Singpurwalla  and  Talbott  found  that 
the  logistics  performance  of  the  C-141  fleet  was  not  improv¬ 
ed;  as  a  result  of  RCM  they  found  some  deterioration.  How¬ 
ever,  they  had  some  criticism  of  their  study  because  they 
did  not  take  some  operational  variables  into  consideration, 
such  as  flying  hours.  The  study  also  ignored  the  fact  that 
the  logistics  performance  variable,  alert  availability  was 
"messy"  (16:2).  The  article  reviewed  here  was  instigated  in 
order  to  rectify  these  deficiencies,  and  to  determine  if 
the  expected  benefits  of  RCM  were  an  improvement  over  the 
previous  system  used  by  the  Air  Force.  The  expected 
benefits  included  reduced  scheduled  maintenance  by 
extending  the  maintenance  intervals,  and  therefore  an 
increase  in  alert  availability. 

The  Box-Jenkins  method  was  used  for  the  transfer  func¬ 
tion  portion  of  the  study.  The  variables,  labeled  X  and  Y, 
represented  the  flying  hours  of  the  fleet  of  C-141 's  and  the 
alert  availability  of  that  fleet  of  aircraft,  respectively. 
The  series  X  was  first  reduced  to  a  white  noise  model  and 
then  because  of  the  "messy"  data  for  the  variable  Y,  it  was 


reduced  to  white  noise  using  an  intervention  analysis  model. 
The  intervention  was  de-fined  as  the  point  in  time  when  the 
RCM  was  introduced  into  the  system,  which  caused  a  change  in 


the  response  of  Y  to  the  input  sequence  (16:4).  The  process 
is  known  as  prewhitening  the  respective  time  series,  which 
is  necessary  be-fore  a  transfer  function  model  can  be  achiev¬ 
ed. 

Singpurwalla  and  Talbott  developed  an  equation  for  their 
final  transfer  function  model  after  they  had  estimated  the 
coefficients  for  the  variables.  The  “...nonzero  coefficients 
associated  with  the  X«’s  imply  that  flying  hours  do  have  an 
effect  on  the  alert  availability  in  a  manner  specified  by 
the  functional  form  of  their  equation  (16:12).  They  found 
that  there  was  a  reduction  in  the  variability  of  the  flying 
hours  as  of  the  inception  of  RCM.  A  concluding  comment  in 
the  article  supports  conclusions  in  some  of  the  articles 
reviewed  earlier.  That  is,  they  pointed  out  that  their  final 
equation  "...supports  the  adage  that,  'the  more  you  fly,  the 
less  you  fail,'  within  limits  (16:12)."  Another  conclusion 
was  that  the  final  model  supported  the  premise  that  RCM  has 
a  tendency  to  reduce  the  alert  availability. 

Use  of  Time  Series  Analysis 

Throughout  this  chapter  comparisons  have  been  made  to 
the  fact  that  there  is  a  "linearity"  issue  and  that  this 
assumption  has  caused  support  problems.  Also  noted  is  the 
fact  that  demands  per  flying  hour  data  is  collected  over 


time  thereby  making  time  a  factor,  which  has  been  ignored. 
Because  time  is  a  factor  there  are  interrelationships  be¬ 
tween  the  two  variables,  demand  and  flying  hours.  This  rela¬ 
tionship  is  not  necessarily  linear,  but  time  series  analysis 
methods  are  the  only  way  to  find  out  what  the  relationship 
really  is.  Once  this  interrelationship  is  found,  a  forecast, 
(which  doesn’t  have  to  assume  a  linear  relationship),  can  be 


CHAPTER  3 


METHODOLOGY 

This  chapter  will  begin  with  an  overview  of  time  series 
analysis,  -Followed  by  detailed  descriptions  o-F  the  implemen¬ 
tation  o-F  these  methods  and  the  data  to  be  analyzed.  A-Fter 
the  overview  the  methods  used  to  prewhiten  the  individual 
series  will  be  described,  then  the  method  -For  combining  the 
prewhitened  models  into  transfer  function  models  will  be 
described.  From  the  transfer  function  model  a  forecasting 
model  can  be  made.  Finally,  a  transfer  function  model  of  one 
pair  of  series,  using  the  original  data  without  prewhiten¬ 
ing,  will  be  shown  to  demonstrate  the  difference  time  series 
methods  make. 

The  data  used  in  this  research  was  provided  by  the  data 
analysis  branch  of  the  F-16  SPO,  (ASD/YPDF) .  This  data  is 
collected  via  a  contract  with  Dynamics  Research  Corporation, 
through  a  central  data  bank  which  is  fed  from  all  F-16  bases 
in  the  Air  Force. 


This  section  will  give  an  overview  of  the  Box-Jenkins 
method  for  the  analysis  of  the  data.  The  method  consists  of 
several  steps,  the  first  of  which  is  using  historical  obser¬ 
vations  of  a  time  series  to  identify  a  tentative  model.  The 


second  step  involves  estimating  the  unknown  parameters  to 
the  tentatively  identified  model.  The  third  step  tests  the 
adequacy  of  the  tentative  model  and  allows  identification  of 
steps  to  improve  it  if  necessary.  When  these  steps  are  fin¬ 
ished,  the  resultant  model  can  be  used  for  forecasting. 

Identification  of  the  tentative  model,  using  the  histor — 
ical  observations,  involves  the  concepts  of  autocorrelation, 
partial  autocorrelation,  cross  correlations  between  the  ser — 
ies  (for  multivariate  series),  and  finding  out  if  the  time 
series  are  stationary  or  non-stat ionary.  Autocorrelation 
describes  the  association  or  mutual  dependence  between  val¬ 
ues  of  the  same  variable  but  at  different  time  periods.  Par¬ 
tial  autocorrelations  are  analogous  to  autocorrelations  in 
that  they  indicate  the  relationship  of  the  values  of  a  time 
series  to  various  time  lagged  values  of  the  same  series. 
However,  they  differ  from  autocorrelation  in  that  they  are 
computed  for  each  time  lag  after  removing  the  effect  of  all 
other  time  lags  on  the  given  time  lag  and  on  the  original 
series.  In  essence,  they  show  the  relative  strength  of  the 
relationship  that  exists  for  varying  time  lags  (13:  247). 
Cross  correlation  is  for  the  multivariate  time  series  what 
the  simple  autocorrelation  is  for  the  univariate  series. 
They  are  standardized  measures,  between  -1  and  +1,  of  the 
association  between  the  present  values  of  a  given  variable 
and  past,  present,  and  future  values  of  another  time  series 
variable.  The  cross  correlations  are  used  to  identify  the 
form  of  the  multivariate  model  (13:587). 


The  residuals  of  the  time  series  must  also  be  checked 


-for  stationarity.  The  patterns  that  the  auto  and  partial 
correlations  produce,  along  with  the  pattern  o-f  the 
residuals,  indicate  the  tentative  model  to  be  used  -for 
estimation  o-f  the  parameters.  This  model  has  the  form  of: 

X*  =  +...+0„X*-p  +  e*  -  +  0,e*-a  - 

••••"  d  q*t-^  (eq.  1) 

if  it  is  stationary.  The  autoregressi ve  'portion  of  the  model 
is  represented  by  the  0i’s,  the  moving  average  portion  is 
represented  by  the  9 t's  and,  the  e4 ’s  represent  the  error 
for  each  observation.  The  subscript  p  and  q  represent  the 
order  of  the  model  for  autoregressi ve  (AR)  and  moving  aver — 
age  (MA)  portions  of  the  model,  respectively.  This  is  con¬ 
ventional  Box— Jenkins  notation.  An  example  of  the  above 
equation  is: 

x*  “  -  0»e*_»  +  e«  <eq.  2) 

where  e  denotes  the  error. 

The  stationarity  can  be  identified  by  observing  the  plot 
of  the  residuals  of  the  model,  if  the  residuals  fluctuate, 
centered  around  the  mean,  it  is  stationary.  If  there  is  not 
a  constant  mean  then  the  model  is  non-st at ionary  and  must  be 
differenced  to  obtain  values  that  can  be  described  by  a  sta¬ 
tionary  process.  The  first  difference  takes  the  form: 

for  t  ■  2, . ,n  (eq.  3) 

Generally,  the  first  difference  of  a  series  is  sufficient  to 
accomplish  a  stationary  condition.  The  second  difference  of 
a  series  takes  the  form: 
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z*  ■  -  Y*  -  2Y*-i  +  Y*-a  for  t  *  3 . ,n 

(eq.  4) 

Occasionally,  it  is  nscsssary  to  take  the  natural  logarithm 
of  a  series  prior  to  differencing  it,  in  order  to  get  it  to 
a  stationary  condition.  The  computer  packages  programmed  for 
time  series  analysis  will  accomplish  all  of  these  calcula¬ 
tions. 

To  identify  the  tentative  model,  after  making  it  sta¬ 
tionary,  one  has  to  look  at  the  autocorrelation  function 
(ACF)  pattern  in  conjunction  with  the  partial  autocorrela¬ 
tion  (PACF)  pattern.  These  patterns  will  indicate  an  autore¬ 
gressive  <AR)  or  a  moving  average  (MA)  process  is  present. 

It  may  be  that  both  are  affecting  the  model,  which  would 
became  an  autor egressi ve  integrated  moving  average  (ARIMA) 
model.  The  integrated  portion  (I)  is  the  level  of  differenc¬ 
ing  required  to  make  the  series  stationary,  it  is  denoted  by 
ad.  The  order  of  p  and  q  is  identified  by  the  pattern  plot¬ 
ted  in  the  PACF.  In  general,  an  AR  model  will  have  an  expo¬ 
nentially  decreasing  pattern  or  an  decreasing  sine  wave  pat¬ 
tern  in  the  ACF  plot.  The  PACF  will  show  a  spike  for  each 
order  of  the  model.  On  the  other  hand,  an  MA  model  will  show 
a  spike  and  then  the  pattern  will  drop  off  sharply  in  the 
ACF.  The  number  of  spikes  will  approximate  the  order  of  the 
model.  The  PACF  for  the  MA  will  show  an  exponentially  de¬ 
creasing  pattern  or  a  decreasing  sine  wave  pattern.  The  or¬ 
der  of  the  model,  <p,d,q),  is  given  the  number  which  corres¬ 


ponds  to  the  number  of  spikes  identified  in  the  ACF  and  PACF 
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and  the  level  of  difference  used 


Once  the  tentative  model  is  identified,  the  values  of 

the  coefficients,  and  must  be  estimated.  This  can  be 

done  by  hand  but  is  very  difficult  for  a  large  time  series. 
The  computer  is  programmed  to  do  this  estimation  and  will  go 
through  an  iterative  process  until  it  gets  the  value  of  the 
mean  squared  error  (MSE)  cannot  be  decreased  further.  At 

this  time  the  final  estimate  for  the  model  is  specified  us¬ 

ing  the  parameters  found  for  d>i.  and  0t.  For  example,  using 
the  general  equation  for  an  ARIMA  (1,0,1)  model  (see  eq.  2), 
if  the  coefficients  are  .30  and  .45  respectively,  the  model 
becomes: 

—  . 45e*_j  +  et  (eq.  5) 

The  test  of  the  adequacy  of  the  model  is  the  next  step. 
This  test  is  simply  checking  to  see  that  residual  differen¬ 
ces  between  the  time  series  and  those  estimated  by  the  model 
are  white  noise.  The  errors  are  represented  by  the  equation: 

(eq.  6) 

The  autocorrelation  coefficient  of  the  residuals,  e*,  can 
then  be  found.  This  will  show  up  as  a  pattern  in  the  residu¬ 
als,  which  is  then  interpreted.  If  none  of  the  autocorrela¬ 
tions  are  significantly  different  than  zero,  the  errors  are 
assumed  to  be  random,  the  residuals  are  white  noise,  and  the 
model  is  adequate.  If  the  errors  are  not  random  then  the 
model  is  still  not  adequate  and  the  above  process,  beginning 
with  the  identification  of  a  tentative  model,  must  be  re¬ 
peated  until  the  final  model  is  all  white  noise  (11:251). 
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All  -Forecasting  methods  assume  some  pattern  or  rela¬ 
tionship  exists  that  can  be  identified  and  used  as  the  ba¬ 
sis  -for  preparing  the  -forecast.  In  order  to  identify  that 
pattern,  the  data  must  be  analyzed,  beginning  with  the  his¬ 
torical  observations  as  they  were  recorded  over  time.  The 
actual  flying  hours  for  each  month  are  assigned  the  indepen¬ 
dent  variable  "X"  and  the  dependent  variable,  HY",  repre¬ 
sents  the  demands  or  removals  per  month.  The  variable  X* , 
Xt-i,  would  represent  the  observed  value  of  actual  flying 
hours  at  the  last  or  latest  observation  and  the  observation 
just  previous  to  that.  The  data  used  in  this  thesis  has  47 
observations  over  a  47  month  time  frame  beginning  with  June 
1979  inclusive  through  April  1983,  so  the  variable  X«._.*7 
would  represent  the  observed  value  for  June  1979.  A  variable 
X«*.w  represents  a  forecast  five  periods  or  months  ahead. 

With  a  transfer  function  the  dependent  series,  Y* ,  can 
be  represented  by  a  linear  operation  on  the  independent  se¬ 
ries,  X%,  plus  a  noise  component  nt  when  Y*,  X*,  and  n*  are 
all  stationary  (18:13).  The  X*  and  Y*  series,  used  in  the 
transfer  function,  are  prewhitened  from  the  original  series. 


The  data  was  analyzed  with  the  TIMES  software  program 
(20).  This  program  is  resident  on  both  the  Cyber  and  the 
Harris/300  computers  at  AFIT. 


Before  the  transfer  function  model  can  be  used,  each 
individual  series  must  be  prewhitened.  An  identity  function 


run  was  made  -For  each  univariate  series  so  that  a  starting 
tentative  model  could  be  identified.  *  The  series  of  actual 
flying  hours  provides  an  illustration  of  the  method  to  iden¬ 
tify  and  prewhiten  a  univariate  series  in  preparation  for 
inclusion  in  the  transfer  function  model.  Figure  1  is  the 
observed  series  for  the  actual  flying  hoursf  it  shows  that 
the  data  is  not  stationary  because  of  the  unequal  deviations 
about  the  mean.  This  non-stationarity  must  be  removed  before 
a  model  can  be  applied.  Figure  2  illustrates  the  scatterplot 
of  the  residuals  after  a  first  order  difference  has  been 
done.  The  model  is  now  stationary,  so  the  ACF  and  PACF  will 
provide  a  proper  graph  to  use  in  determining  a  tentative 
model.  The  autocorrelation  function  (ACF),  Figure  3  ,  and 
the  partial  autocorrelation  function  (PACF),  Figure  4,  are 
examined  to  determine  the  order  of  the  model.  An  autoregres¬ 
sive  integrated  moving  average  (AR1MA)  model  is  specified 
which  will  remove  all  the  correlation  within  the  variables 
and  leave  only  random  noise  for  the  prewhitened  models.  The 
residuals  from  this  prewhitened  series  provide  the  input  to 
the  transfer  function  model. 
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Scatterplot  o-f  the  unaltered  scries  (actual  f 
juts).  Tha  zero  lina  with  tha  MMM*s  reprasents  the  me 
>  series.  Tha  plot  shows  how  tha  series  deviates 
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Figure  2 

Deviations  from  the  Mean  -  Stationary 


Scatterplot  of  the  residuals  after  first  order  differen 
cing.  The  residuals  are  now  stationary  about  the  mean. 
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Figure  3 

Autocorrelation  Function 

The  negative  spike  at  lag  1  is  significant  because  it  is 
greater  than  two  standard  deviations  (+/-. 300). 
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Figure  4 

Partial  Autocorrelation  Function 


The  autocorrelation  function  and  this  partial  autocor¬ 
relation  function  have  significant  spikes  which  are  used  to 
determine  the  tentative  model. 
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The  original  graphs  of  the  ACF  and  the  PACF  indicated 
that  the  tentative  model  should  be  an  auto  regressive  order 
one.  The  patterns  o-f  the  series  after  differencing  (charts  3 
and  4)  have  significant  spikes  at  lag  1  in  the  ACF  and  lags 
1  and  2  in  the  PACF.  A  significant  spike  is  one  that  extends 
beyond  two  standard  deviations  (in  this  case  beyond  ■*■.300 
and  —.300).  The  two  spikes  in  the  PACF  are  indicative  of  an 
order  2  model.  The  pattern  in  the  ACF  also  shows  a  dampening 
sinusoidal  pattern  further  indicating  an  AR  function.  The 
process  of  identifying  a  final  model  is  an  iterative  one,  in 
this  case  an  ARIMA  (2,1,0),  as  is  indicated  in  the  printout 
of  the  identification  run,  is  the  most  probable  model  to 
start  with.  An  ARIMA  (2,1,0)  model  is  an  AR  order  2  with  a 
difference  of  order  1  and  no  MA  function. 

When  the  correct,  most  parsimonious,  model  is  found  the 
residual  mean  square  will  be  at  the  smallest  possible  value. 
Also  the  Chi  Square  value  shown  in  Figure  3  will  be  used  to 
verify  that  the  model  is  white  noise.  The  residual  mean 
square  for  the  identification  model  with  the  first  differ¬ 
ence  taken  to  make  the  series  stationary  is  found  in  Figure 
5.  It  should  also  be  noted  that  the  data  was  transformed 
with  a  logarithmic  function.  This  was  done  because  the  data 
for  the  actual  flying  hours  are  very  large  numbers,  compared 
to  the  demand  data. 

Figure  b  is  a  power  spectrum  of  the  actual  flying  hour 
residuals  after  taking  a  first  difference.  From  this  graph. 


seasonal ity  and  the  fit  of  the  data  can  be  determined. 
Ideally,  the  power  spectrum  of  the  residuals  should  be  as 
flat  as  possible,  that  is,  a  horizontal  line  drawn  through 
the  spectrum  should  be  centered  on  the  SSS’s,  thereby  indi¬ 
cating  that  the  residuals  are  white  noise. 


Figure  5 

Summary  of  Model 
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The  Final  Prewhitened  Univariate  Model 

The  -final  model  for  the  actual  -flying  hours  is  an  ARIMA 
(2,1,2).  This  model  can  be  represented  by  the  equation: 

(l-0iB-fc,B*)(l-B)Z*  -  (1-&B— 02B=>a*  (eq.  7) 

The  backshi-ft  operator,  B,  is  a  Box-Jenkins  notation 
that  is  used  to  simplify  the  writing  o-f  an  equation.  In  the 
above  equation  the  backshi-ft  operator  Bs  is  the  equivalent 
o-f  Z*_»  and  the  random  error  is  represented  by  a*.  The  d4h 
difference  of  Z*  is  represented  by  (l-BJ^Z*. 

Figure  7  is  a  periodogram  of  the  ARIMA  (2,1,2)  model. 
This  represents  the  "fit"  of  the  series,  a  perfect  fit  would 
follow  a  diagonal  line  through  the  center  of  the  graph.  This 
fit  is  as  close  as  possible  with  the  number  of  observations 
available.  It  should  be  noted  here  that  as  more  observations 
are  accumulated  it  may  be  possible  to  improve  this  model  of 
the  actual  flying  hours  time  series. 

Figures  8  and  9  are  the  ACF  and  PACF  functions  of  the 
ARIMA  (2,1,2)  model.  The  spikes  that  were  present  in  the 
identification  model  are  now  gone  and  all  lags  are  within 
the  two  standard  deviations.  Note  that  one  standard  error 
can  be  found  in  Figure  8,  it  is  printed  out  for  each  lag  and 
listed  in  the  row  labeled  ST.E.  Two  standard  deviations  are 
needed  to  determine  if  the  pattern  is  white  noise,  these  are 
found  by  either  multiplying  the  figure  given  in  the  ST.  E. 
or  by  Bartletc  s  approximation.  Bartlett's  approximation  for 
two  standard  errors  is  equal  to  2  times  the  square  root  of 
the  reciprocal  of  the  number  of  observations  (17). 
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The  Chi -Squared  test  is  used  to  test  whether  the  mean 
•values  o-f  the  autocorrelations  is  sign if icantly  different 
than  zero.  If  the  value  of  the  Chi-Square  test  is  smaller 
than  the  corresponding  value  from  a  table  of  computed  Chi- 
Square  values,  then  the  data  are  random,  without  any  pattern 
(13:583).  Figure  8  gives  a  computed  Chi -Squared  value  of 
6.0612  with  11  degrees  of  freedom,  this  is  compared  with  a 
tabled  value  of  17.2750  at  X*.  Thus  the  data  are  not 
autocorr el ated . 


•f.3N  ♦ 


Figure  10  is  a  sunnary  of  the  final  nodel  for  the  actual 
flying  hours.  This  nodel  will  be  used  in  the  transfer  func¬ 
tion  nodel s  where  it  is  combined  with  the  other  univariate 
models  to  find  the  cross  correlation,  or  "between"  correla¬ 
tion,  so  that  a  forecast  can  be  accomplished.  Note  that  the 
data  was  transformed  with  a  logarithmic  function.  The  pur — 
pose  of  this  is  to  reduce  the  magnitude  of  the  monthly  ob¬ 
servations,  thus,  making  the  residuals  smaller  in  their  mag¬ 
nitude.  This  must  be  done  in  order  to  find  the  impulse  re¬ 
sponse  weights  during  the  identity  of  the  transfer  function. 
The  impulse  response  is  a  series  of  weights  in  a  multivari¬ 
ate  (ARIMA)  model  relating  the  independent,  input,  series  to 
the  dependent,  output,  series.  The  value  of  each  weight  de¬ 
termines  the  relative  importance  of  each  residual  to  the 
others,  a  large  weight  will  cause  an  impulse  in  the  pattern 
of  the  residuals.  This  pattern,  and  the  weights  are  useful 1 
in  identifying  the  model.  The  following  section  will  give 
more  detail  on  this  process. 

The  first  part  of  Figure  10  is  where  the  information 
about  the  coefficients  of  the  AR  and  HA  parameters  are 
given.  By  inserting  this  information  into  equation  7,  the 
final  model  for  the  actual  flying  hours  becomes: 

(1  +  1.0122B  -  . 0084937B® Ml  -  BJZ*  - 

(1  +  . 78672B  -. 43365Ba) a« 

The  residual  mean  square  for  this  model  is  .0086986, 
which  was  the  smallest  found  in  the  iterative  process  in 


developing  the  model.  Finally,  the  power  spectrum  o-f  the 
residuals  is  indicative  o-f  white  noise.  As  noted  earlier, 
this  spectrum  ideally  should  have  even  distribution  o-f  the 
SSS*s  about  a  horizontal  line  drawn  through  the  graph.  The 
distribution  in  this  case  is  satisfactory  and  represents  the 
best  possible  fit  of  the  model  to  the  data.  A  line  drawn  at 
approximately  the  -2.000  level  will  have  an  equal  distribu¬ 
tion  above  and  below  it. 

Appendix  B  contains  summary  printouts  for  the  balance  of 
the  univariate  models.  These  are  arranged  in  the  same  order 
as  the  work  unit  codes  (WUC’s)  in  Appendix  A.  Equations  for 
each  of  the  five  WUC’s  and  the  corresponding  ARIMA  model 

are: 

WUC:42GB0  ARIMA  (0,1,1) 

<1  -  0X1  -  B)Z*  -  (1  -  . S3800B ) a.*. 

WUCM4FB0  ARIMA  (2,1,1) 

(1  +  .0068762B  +  . 47673B*) (1  -  B)Z«  =  (1  -  ,75544B)a*. 

WUC: 24EBA  ARIMA  (2,1,1) 

(1  +. 12117B  +  .31373b*) (1  -  B)Z*  »  (1  -  .52364B)a« 

WUC: 46AF0  ARIMA  (2,1,1) 

(1  +  . 10442B  +  . 38446B*) (1  -  B)Z*  ■  (1  -  .76877B)a« 

WUC: 46AN0  ARIMA  (2,1,0) 

(1  +  . 70241B  +  . 281 00B*) ( 1  -  B) Z*  »  (1  -  0)a* 

The  fit  of  the  individual  univariate  models  can  be  seen  in 

the  power  spectrum  for  each  model. 
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Figure  10 

Summary  of  Actual  Flying  Hours  Model  -  ARIMA  (2,1,2) 
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The  Multivariate  Model 

The  -first  step,  after  the  prewhitening  of  each  univari¬ 
ate  series,  is  to  identify  the  class  of  model  for  the  trans¬ 
fer  function.  The  pairs  of  time  series  models  have  been  dif¬ 
ferenced  and  made  stationary.  The  output  series,  Y*  (the 
demands)  and  the  input  series,  X*  (the  actual  flying  hours), 
are  now  represented  by  the  equation; 

Y*  *  (v«+ViBx+vaBa+. . . ) X«+e«  (eq.  8) 

where  Y *  is  described  as  the  weighted  sum  of  lagged  values 
from  both  the  input  and  output  series  plus  random  error  (13: 
379).  The  impulse  response  weights,  v^s,  are  difficult  to 
estimate  and  therefore,  the  above  equation  is  not  an  effi¬ 
cient  parameteri2ation  (15:380). 

According  to  Box  and  Jenkins,  the  model  formula¬ 
tion  lacks  "parsimony."  The  number  of  parameter  esti¬ 
mates  could  result  in  an  unmtmblm  estimation  problem 
CIS: 3791. 

A  more  general  and  parsimonious  equation  for  the  trans¬ 
fer  function  is: 

(1-iiB-.  .  .-drB«-)Y*  -  <W«-OK B- - -W.B-)  X,._*+Nt.  (eq.  9) 

where  r  is  equivalent  to  the  number  of  cogent  past  Y«  val¬ 
ues,  s  is  equivalent  to  the  number  of  past  X*  values,  b  is 
the  lag  factor  (the  number  of  periods  until  the  effect  of  a 
change  in  X*  affects  Y*),  and  N*  represents  the  noise  fac¬ 
tor.  The  Y«  is  the  same  as  found  in  eq.  8.  This  model  is 
referred  to  as  an  (r,s,b)  model. 

The  impulse  response  weights,  v<’s,  are  used  to  provide 
guidelines  to  the  identification  of  the  appropriate  values 
of  the  (r,s,b)  model.  Box  and  Jenkins  demonstrate  that  the 


impulse  response  weight  v*  is  simply  lag  k  cross-correlation 
of  the  a *  and  series,  <k),  multiplied  by  a  scale  fac¬ 
tor  (6:  380).  The  cc*  and  f3*  series  are  the  prewhitened  input 
and  output  series,  respecti vely.  The  equation  for  finding  v* 
is: 


v*  =  Pafi< k)  k=0, 1,2, _ _  (eq.  10) 


Equation  (10)  suggests  that  the  series  referred  to 
as  "prewhitened"  retain  information  on  (Y*,  X*)  inter¬ 
relationships  without  the  sometimes  confounding  infor — 
mation  on  intrarelationships  in  the  individual  series. 
Prewhitening  removes  the  latter  information  to  illu¬ 
minate  the  former  CIS: 3823. 

The  impulse  response  weights,  vk,  will  have  groupings 
which  give  clues  to  identify  the  correct  values  for  (r,s,b). 
The  quotations  below  relate  to  eq.  8  &  9  above. 

The  clue  for  b  is: 

for  k  ■  0,1,..., b-1,  the  right  side  coefficient  of  B* 
is  0.  Thus,  v«»v1*!...»vb-,«0. 

The  clue  for  r  is: 

for  k>b+s,  the  right  side  coefficient  of  B*  is  0.  Ex¬ 
pansion  of  the  left  side  then  implies  that  vK-iavk_t- 

(JrVK-r-*0.  vh  generated  from  an  equation  of  this 
type  would  illustrate  patterned  behavior.  For  r»l,  a 
pattern  of  exponential  decrease  would  be  seen  (with  v„ 
alternating  in  sign  if  ji<0).  For  r«2,  a  typical  pat¬ 
tern  would  be  that  of  a  damped  sine  wave,  although 
exponential  decrease  is  also  likely.  Thus,  the  value 
of  r  can  be  identified  through  association  of  impulse 
response  weight  patterns  that  are  exactly  equivalent 
to  the  autocorrelation  patterns  associated  with  the 
same  value  of  the  autoregressi ve  order  p  in  a  univari¬ 
ate  time  series  model.  An  additional  aid  in  identify¬ 
ing  r  may  exist  in  the  visible  grouping  of  r  "initial 
values"  or  "startup  values"  for  the  equation  given 
above.  These  values  will  be  the  v*  that  appear  in  the 
•xpression  for  Vb...t,  i.e.,  Vb+.+i , . .  . ,Vb-.«-t ,  Vu... 

Often,  however,  these  r  startup  values  will  seem  to  be 
part  of  the  pattern  in  the  v„  for  k>b+s. 

The  clue  for  s  is: 

the  pattern  referred  to  in  the  clue  for  the  identifi- 


cation  of  r  begins  with  the  impulse  response  weight 
Vt3 i .  Establishing  the  beginning  of  this  pattern 
establishes  the  value  of  b+s+1.  Given  previous  identi¬ 
fication  of  b,  one  may  deduce  a  most  likely  value  for 


The  Transfer  Function  Model 

The  transfer  function  model  to  be  used  for  this  example 
is  of  the  WUC: 14FB0.  Appendix  C  contains  the  remainder  of 
the  transfer  function  models. 

The  identity  run  for  the  model  used  in  this  example 
gives  the  following  data  to  use  in  identifying  a  tentative 
model.  The  effective  number  of  observations  for  cross  cor¬ 
relation  calculation  is  44. 


Table  1 


Model  Information 


Mean  of  Series  1  ■ 
St.  Dev.  of  Series  1  “ 
Mean  of  Series  2  ■ 
St.  Dev.  of  Series  2  » 


0.75860E— 11 
0. 23736E+00 
0. 11038E-10 
0. 84308E+00 


The  mean  values  in  Table  1  are  used  with  eq.  9  to  esti¬ 


mate  the  v*  values  at  each  lag.  This  is  accomplished  by  the 


TIMES  program  which  saves  the  effort  of  calculating  manual¬ 
ly.  Table  2  was  extracted  from  the  identity  run  to  show  the 


calculated  impulse  response  weights,  v*.  Note  that  lag  6  is 
a  large  negative  value,  this  shows  up  in  the  plot  in  Figure 
11. 
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Table  2. 


Estimated  Impulse  Response  Weights  v(k>. 

!<  V  <K) 

0  0.777 

1  -0.075 

2  0. 186 

3  0.110 

4  -0.035 

5  0.501 

6  -1.391 

7  0.617 

8  0.416 


Table  3  is  the  cross  correlation  values  that  are  plotted 
in  Figure  11.  Series  1  data  is  plotted  on  the  right  side  of 
the  graph. 


Table  3 


Cross  Correlations 

Series  1  -  Prewhitened  Actual  Flying  Hours 
Series  2  -  Prewhitened  WUC: 14FB0  El  Comp  Assy 


Number  of  Lags 
on  Series  1 


Cross  Number  of  Lags 

Correlations  on  Series  2 


Cross 

Correlations 


0 

1 

2 

3 

4 

5 

6 

7 

8 


0.237  0 

-0.023  1 

0. 057  2 

0.033  3 

-0.011  4 

0. 153  5 

-0. 425  6 

0. 188  7 

0. 127  8 


0.237 
-0. 128 
-0.083 
0.073 
0.254 
-0. 180 
-0.046 
0.081 
0.044 


Figure  11  represents  the  cross  correlations  between  the 
pair  of  time  series.  The  right  side  of  this  plot,  (right  of 
zero),  is  important  in  that  it  is  the  side  that  is  used  to 
identify  the  final  transfer  function  model.  When  the  model 


is  correct  this  side  will  have  all  of  its  spikes  within  two 
standard  errors. 


Figure  11 

Cross  Correlation  Plot 
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Figure  12  is  the  estimated  impulse  response  -for  this 


model.  This  -figure  is  a  plot  o-f  the  v*  weights  of  the  esti¬ 
mated  cross  correlation.  With  this  graphical  portrayal  and 
Figure  11  one  can  estimate  the  values  for  r,  s,  and  b.  It  is 
obvious  that  there  is  not  a  time  lag  between  the  input  and 
output  series  because  of  the  large  spike  at  lag  0  indicating 
that  V(g»  is  not  equal  to  zero,  therefore,  b  =  0  for  this  mod¬ 
el.  The  r  value  is  also  zero  because  there  is  no  apparent 
pattern  in  the  cross  correlation  function.  The  s  value  is 
all  that  remains  to  figure  out.  This  is  easily  accomplished 
by  looking  at  the  cross  correlation  produced  by  the  identity 
run  (Figure  11).  The  lags  that  have  spikes  which  are  less 
than  two  standard  errors  long  have  an  expected  value  of  ze¬ 
ro,  therefore,  they  are  not  significantly  correlated.  The 
lags  with  spikes  longer  than  two  standard  deviations  are 
produced  by  significant  inpulse  response  weights  from  which 
the  s  value  can  be  found.  The  standard  error  can  be  found  by 
taking  the  square  root  of  the  reciprocal  of  the  n  observa¬ 
tions  (6:382). 
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Figure  12 

Estimated  Impulse  Weights 


This  figure  represents  the  estimated  impulse  response 
that  are  used  to  determine  the  tentative  transfer  function 
model 
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Figure  13  is  a  summary  o-f  the  transfer  function  model. 
The  same  type  of  information  as  found  in  the  summaries  of 

the  univariate  models  is  found  here.  In  addition,  the  opti¬ 
mum  value  of  b  is  given,  any  noise  model  parameters  are  giv¬ 
en  (in  this  case,  since  both  univariate  models  were  pre— 
whitened,  no  noise  parameters  were  required).  The  estimated 
values  for  the  transfer  function  parameters  are  the  most 
important  information  found  in  this  Figure. 

Figure  13 

Transfer  Function  Summary 

SUMMARY  OF  MODEL 

DATA  -  X  -  ACTUAL  FLY 1 NO  HOURS  44  OBSERVATIONS 

Y  -  MUCI 14FDS  EL  COMP  ASSY 

DIPPERENCXNS  ON  Y  -  NONE 

DIFFERING INS  ON  X  —  NONE 

NOISE  MODEL  PARAMETERS 
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TRANSFER  FUNCTION  PARAMETERS 


1  INPUT  LAS  1  S  ■ 7SSS4E+S#  143SSE+SA  . 17349E+SI 

2  INPUT  LAS  1  2  -.144*3E*S#  -.U494E+SI  .S2137E*## 

3  INPUT  LAS  I  4  . 79739E+SS  -.29939E+SS  .177#*E*S1 

4  INPUT  LAB  1  9  . 4S72SE-S1  ~.9t92SE+S#  . 1S927E+S1 

9  INPUT  LAS  I  A  .  17S23E*#!  .7344SE«SS  .244478*01 

4  INPUT  LAS  I  7  -• 44S43E+SS  -.139398*01  .497#2E*S# 


OPTIMUM  VALUE  OF  D  IS  # 


RESIDUAL  SUM  OF  S8UARES 


.140978*02  31  D.F. 


U: 
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RESIDUAL  STANDARD  ERROR 


From  the  data  in  Figure  13,  an  equation  for  the  transfer 
function  model  can  be  written  by  inserting  the  numbers  into 
eq.  9.  The  equation  for  this  function  is: 

Yt  =  (.79554  +  . 16430B3  -  .757358*  -  .068728Ba  -  .0170238*  + 
.  44843B7’)  Xt  +  Nt . 

The  periodogram  in  Figure  14  represents  the  fit  of  the 
transfer  function  model  within  .1  probability  limits. 


Figure  14 

Transfer  Function  Periodogram 
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The  histogram  in  Figure  IS  represents  the  distribution 
of  the  residuals  in  the  transfer  function  model.  This  histo¬ 
gram  is  a  near  perfect  normal  distribution,  the  histograms 
for  the  balance  of  the  transfer  function  models  are  found  in 
Appendix  C.  The  others  transfer  function  models  do  not  have 
as  perfect  a  distribution  as  this  one  but  they  are  still 
well  within  the  expected  range  of  distribution. 


Figure  15 


Histogram  for  WUC:  14FB0  Transfer  Function. 
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The  cross  correlation  in  Figure  16  shows  that  none  o-f 
the  lags  on  the  right  side  have  spikes  significantly 
different  than  zero.  This  model  has  all  of  the  between 
structure  correlation  removed  and  is  now  ready  to  be  used  in 
forecasting.  A  review  of  all  the  cross  correlation  figures 
in  Appendix  C  will  show  that  all  of  the  final  transfer 
function  models  have  only  white  noise  on  the  right  side  of 
their  cross  correlation  plots.  The  equation  for  Figure  16 
has  been  written  above,  equations  for  the  balance  of  the 
transfer  function  models,  by  WUC  and  type  of  model  are: 

WUC:42GB0  (0,8,0) 

Y*  =  (-4.8069  -  3.5346B  +  4.9247B*  -  5.2819B^  +  1.9522B°)X* 

+  N*. 

WUC: 24EBA  (2,4,0) 

(1  -  • 077362B  +  . 13579B2) Y*  =  (.79689  -  1.4642B  +  .8B371B=  - 

2.6610B^)X*  +  N-t . 

WUC: 46AF0  (1,6,0) 

(1  -  . 35718B) Y*  =  (2.3915B  -  1.3710B0  +  1.3305B*>)X*  +  N*. 

(1,5,0) 

=  (.52943  -  1.3678B*  +  .61573B*  - 
2.  2439B55)  Xt  +  N*. 


WUC: 46AN0 

(1  +  .  45490B)  Y-*-. 


Figure  16 

Cross  Correlations  -for  WUC:  14FB0 
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i«f«r  Function  Using  Raw  Data 
The  present  method  of  forecasting  requirements  is  to  use 
the  linear  function,  Y*  «  a«,  +  axXt  +  et.  When  this  type  of 
equation  is  used,  the  time  factor  is  ignored,  the  possibili¬ 
ty  of  seasonal  factors  is  ignored,  and  it’s  assumed  that  the 
model  is  adequate.  This  section  will  show  by  comparison,  the 
difference  in  using  raw  data,  assumed  linear,  and  the  time 
series  method  of  the  previous  sections.  To  do  this,  the  data 
was  input  to  the  TIMES  program  as  a  transfer  function  which 
compares  two  series  of  data  to  determine  the  between  series 
correlations.  In  the  section  on  univariate  time  series  it 
was  pointed  out  that  the  data  was  non-st at ionary  and  there¬ 
fore  had  to  be  differenced  prior  to  finding  a  correct  model. 
The  same  set  of  time  series  used  in  the  multivariate  section 
was  input,  in  it’s  raw  form,  for  this  section.  This  is  the 
form  that  would  be  used  in  the  present  forecasting  methods, 
i.e. ,  just  as  it  was  collected.  Compare  the  figures  in  this 
section  to  those  in  the  previous  section  to  see  that  the 
methods  of  time  series  analysis  do  a  better  job  of  fitting 
the  data.  Figure  17  is  a  plot  of  the  residuals  about  the 
mean  (MMM) ,  this  plot  comes  from  the  run  which  will  be  call¬ 
ed  the  ’old  method’  in  the  figures  and  text.  Note  that  it  is 
not  stationary,  that  is,  it  does  not  vary  uniformly  about 
the  mean,  but  trends  across  the  mean.  Compare  this  plot  with 
Figure  18,  which  is  from  the  run  that  produced  the  figures 
for  the  previous  section. 
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Figure  17 


Deviations  from  the  Mean 
The  ’Old  Method' 
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Figure  18 

Deviations  -from  the  Mean 
Time  Series  Model 
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The  periodogram  below  should  be  compared  to  Figure  14. 
Note  that  the  residuals  do  not  fit  the  graph  at  all  in  this 
periodogram,  demonstrating  that  the  residuals  from  the  ’Old 
Method’  are  still  autocorrel ated 


Figure  19 

Periodogram  for  the  ’Old  Method’ 
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This  histogram  is  compared  to  Figure  15,  this  plot  is 
skewed  to  the  right  and  Figure  15  is  near  perfect.  The  his¬ 
togram  is  not  as  good  of  an  indicator  of  fit  as  the  periodo- 
gram  and  spectrogram  and  should  not  be  relied  on,  exclusive 
of  the  other  indicators. 

Figure  20 

Hi stogram 
’Old  Method’ 


CELL  WIDTH  .37541*02 

3.90  * - - — - 

S  i 


2.8#  ♦/ 


/////✓✓///✓///////////✓/ I 
Ittlttlt/ftttt/ttt/tttttX 
//////////////✓//////✓// X 
/////////////////////✓//I 
/✓//✓////✓/✓//////////✓/ I 
//////////////////////// I 
/////✓///✓//////////////l 

tnttntt+tutntt  /♦///! - ♦- — - ♦ — — - -♦ - - - - — — 

//✓✓✓/✓//////✓✓✓//✓✓//// i 

/✓////////////////////✓/ x 

/////////✓///✓✓///✓/✓/// i 
//✓///////✓/✓✓✓✓✓//✓✓/// i 
//✓////////////////✓✓✓// i 
///✓//✓//////////////✓✓/ i 
//////✓/✓/✓///✓/✓/✓/✓/// i 
//////////✓///////////// i 
/////✓/////✓/✓//////////i 

2.10  ♦✓//✓/////♦/////////♦✓//! - - — — — ♦* — * - - — — - ♦— — — — — 

////////////////////////i 

////////////////////////x 
////////////////////////x 
//ntf/mntftttfttitn  I 

//✓//✓✓///✓///✓//////✓/✓ i 
////////////////////////% 

////////////////////////x 

✓✓✓///////✓///////////// i 

1.40  +/////////  +  /////////  +  ///1 — ♦ - - - - - ♦ - ♦ — — - - - 

//✓///✓///✓///////✓/✓✓// I 
////////////////////////t 
////////////////////////l 
////////////////////////l 
////////////////////////l 

////////////////////////X////////////////////////X 
✓//✓//✓///✓/ /✓/✓///////✓!/✓///////✓/////✓//// ////I 
////f///////t///////////l///////////////////f////X 
//////////////////////// I ✓✓///////////✓/////// ///l 

.70  ♦/////✓///♦/✓✓/✓✓✓✓✓♦/✓/ !///✓/♦/////////♦//// ////!♦ - - - ♦ - - — - 

//////////////////✓/✓/// I ////////✓//////✓//////// I 
//////////////////////// X////////////////////////X 
////////////////////////X////////////////////////1 
////////////////////////X////////////////////////X 
////</✓ ///////✓//////// /!/✓///////////////✓//////! 

//////////////////////// 1////S///////////////////X ///////////////////// ///l 
//////✓///<///✓///✓///// !/✓//////✓///////////////!///////////////////// ///I 
ttt/uttttntntttttttnxnttftitttitttttt/uttnxtt/ttutitttfttnutttttx 
n/tututtnuuttuitt  x//// /✓///// /////////////x//////////////////////// 1 
0.00  ♦✓///////✓♦/✓////✓✓/♦/✓/ I /////♦/////////♦////////!♦///////✓/♦//////// /♦///!- 
-3.000  8  I  -.000  1.300  3.000  4.300  4.000  7.300 


7.000 


The  linear  equation,  Y«  «  a«  +  a*X*  +  e«f  is  put  in  the 
form  of  the  transfer  function  where  Y*  ■  The  factor 

a«  takes  the  expected  value  zero.  The  estimated  value  of 
input  lag  1,  parameter  order  0,  in  Figure  21,  .00014143, 

represents  (Je> .  Compare  also  the  residual  mean  square  value 
in  both  figures.  The  transfer  function  model  portrayed  here 
is  a  (0,1,0).  When  the  time  factor  of  the  data  is  ignored, 
this  is  the  type  of  model  and  fit  that  will  occur. 

Figure  21 

Summary  of  Model 
'□Id  Method' 
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Figure  22  demonstrates  that  there  is  correlation  between 


the  series. 


Compare  this  -figure  to  Figure  16. 


Figure  22 


Cross  Correlations 
'Old  Method* 
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CHAPTER  4 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  last  two  sections  in  Chapter  3  demonstrated  multi¬ 
variate  time  series  analysis  with  the  data  fitting  a  model 
in  the  first  section  and  the  data  left  ’raw’  in  the  last 
section.  The  two  sections  used  the  same  pair  of  time  series 
data,  but  it  can  easily  be  seen  that  when  the  data  is  left 
in  it's  raw  state  that  it  does  not  fit  into  a  model.  This 
state  is  similar  to  that  which  is  presently  used  to  forecast 
spares  requirements.  This  chapter  will  draw  some  conclusions 
about  the  study  and  make  some  recommendations  about  how  it 
can  be  improved  and  incorporated  into  the  Air  Force  system 
of  spares  forecasting. 

Findings 

An  examination  of  the  transfer  function  model,  in  the 
section  titled  The  Transfer  Function  Model,  of  Chapter  3 
shows  that  there  is  a  seven  period  time  lag  for  the  depen¬ 
dent  variable,  WUC:14FB0,  before  the  independent  variable, 
actual  flying  hours,  has  an  effect  on  it.  It  was  also  shown 
that  both  variables  were  non-st at ionary  and  had  to  be  made 
stationary  before  they  could  be  used.  Both  series  had  auto¬ 
regressive  and  moving  average  components  and  the  interrela- 


tionships  within  each  series  had  to  be  dealt  with  before  the 
intrarelationships  between  them  could  be  -Found. 

These  -Facts  argue  against  the  linearity  assumption  dis¬ 
cussed  in  Chapters  1  and  2.  Comparing  the  trans-Fer  -Function 
model  with  the  ’Old  Method’  model,  described  in  Chapter  3, 
demonstrates  a  clear  difference  between  the  two  which  would 
produce  very  different  forecasts.  The  ’Old  Method’  model 
does  not  remove  the  correlation  between  and  within  the  se¬ 
ries,  as  demonstrated  in  Figure  19.  This  figure  demonstrates 
that  the  relationship  is  not  a  linear  one,  and  should  not  be 
used  without  first  considering  the  elements  which  produced 
the  much  better  fit  shown  in  Figure  14.  Without  first  remov¬ 
ing  the  "within  and  between  structure"  relationships,  it  is 
not  possible  to  accurately  forecast  with  the  raw  data. 


Time  series  analysis  should  use  at  least  50  and  prefer — 
ably  100  observations  to  realize  the  full  benefit  of  its 
powerful  ability  to  produce  accurate  forecasts.  This  thesis 
used  selected  WUC’s  from  the  F-16  which  did  not  have  that 
many  observations  because  of  it’s  recent  introduction  into 
the  Air  Force  inventory.  As  this  aircraft  become  older  the 
methods  demonstrated  in  this  thesis  will  become  easier  to 
use.  That  is,  with  more  observations  to  work  with,  the  pat¬ 
terns  produced  with  the  TIMES  program  are  more  readily  iden¬ 
tifiable.  This  would  make  the  modeling  process  less  time 
consuming  and  enable  more  precise  identification  of  the  re- 


lationships,  thereby  reducing  the  iterations  required  to 
arrive  at  the  correct  model.*  A  rule  of  thumb  in  any  statis¬ 
tical  effort  says  that  "the  more  observations  available,  the 


more  accurate  the  findings  will  be,"  this  applies  equally 
wel 1  to  time  series  analysis. 

Conclusions 

This  thesis  has  attempted  to  demonstrate  that  using  the 
"linearity  assumption"  in  the  business  of  forecasting  can 
lead  to  errors,  because  this  assumption  ignores  the  time  fac¬ 
tor  inherent  in  the  data  collection  process.  Ignorance  of 
this  factor  produces  forecasting  models  such  as  demonstrated 
in  the  ’Old  Method’  section  of  Chapter  3.  This  model  clearly 
has  cross  correlation  present  between  the  variables  which 
will  cause  errors  in  a  forecast.  The  thrust  of  the  thesis 
effort  has  been  to  show  that  time  series  analysis  will  pro¬ 
duce  models  better  suited  to  forecasting  future  spares  re¬ 
quirements.  When  the  models  produced  with  time  series  meth¬ 
ods  are  compared  to  a  model  of  the  same  data,  produced  with¬ 
out  the  benefit  of  time  series  methods,  then  time  series 
forecasting  becomes  a  more  intuitively  appealing  method  of 
estimating  what  the  future  spares  requirements  might  be. 

With  all  the  time  series  analysis  programs  readily  a- 
vailable,  and  the  increasing  availability  of  computer  re¬ 
sources  to  the  IM’s,  the  assumption  that  demands  per  flying 
hour  is  just  a  linear  ratio  need  no  longer  be  made.  In  addi¬ 
tion,  the  improved  data  collection  systems  being  introduced 


enable  the  IM  to  get  all  the  data  needed  to  make  more  accur — 
ate  forecasts  for  the  spares  requirements.  This  could  lend 
itself  to  dollar  savings  as  well  as  improved  readiness  and 
support  of  weapon  systems. 

The  TIMES  software  package  belongs  to  the  Air  Force  and 
can  be  booted  to  any  main  frame  computer  that  is  capable  of 
compiling  FORTRAN  ver.  5  (20).  There  is  also  a  software  pro¬ 
gram  called  SIBYL/RUNNER  (13)  which  is  available  on  the  Har¬ 
ris  S00.  This  is  a  more  user  friendly  time  series  package 
although  I  do  not  believe  it  is  as  powerful  as  the  TIMES 
program.  The  SPSS  and  BMDP  packages  also  include  time  series 
analysis  programs. 

Once  a  model  is  developed,  the  program  allows  it  to  be 
continually  updated  with  future  observations.  This  implies 
that  once  the  effort  to  develop  models  for  each  WUC  or  LRU 
is  accomplished,  it  is  a  simple  matter  of  adding  the  latest 
observations  to  update  the  forecast.  All  of  the  software 
packages  allow  the  user  to  specify  forecasts  for  future  per¬ 
iods  which  they  define.  Updating  the  model  periodically  will 
assure  that  long  term  forecasts  remain  accurate. 


Recommendat ions 

The  models  in  this  thesis  should  be  verified  with  addit¬ 
ional  observations  as  they  become  available.  Also,  the  com¬ 
puter  resources  available  to  the  Item  Manager  should  have 
the  TIMES  software  package  implemented  so  that  they  can  be¬ 


gin  using  it  to  aid  them  in  their  support  of  the  various 


weapon  systems.  The  initial  effort  to  develop  -forecasting 
models  will  pay  -for  itsel-f  in  future  time  saved,  improved 
support  and  readiness  of  the  individual  systems.  The  alloca¬ 
tion  of  funds  for  spares  is  never  enough,  especially  if  the 
spares  purchased  are  based  on  an  assumption  that  their  util¬ 
ization  is  based  on  a  straight  linear  relationship  with  act¬ 
ual  flying  hours.  There  is  a  lot  of  concern  that  support 
based  on  this  assumption  in  peacetime  will  be  short  of  some 
requirements  in  a  wartime  situation.  This  concern  is  great¬ 
est  where  the  more  expensive  LRU’s  are  the  required  spare 
(7$8$9$ 10$ 11$ 12) .  The  required  spare  should  be  available 
when  needed,  regardless  of  its  cost,  and  using  time  series 
methods  to  forecast  should  increase  the  probability  that  it 
will  be  there,  without  increasing  the  entire  inventory. 


Figure  A1 


The  Data  Used  In  Time  Series  Models 


P 


V 

y 

j 

j 

!t 


Actual 

Flying  Removals  -for  Work  Unit  Codes 


Date 

Hours 

42GB0 

14FB0 

24EBA 

46AF0 

46AN0 

7906 

411.55 

0 

2 

2 

0 

0 

7907 

552. 35 

0 

2 

0 

0 

0 

7908 

516.50 

0 

5 

0 

0 

0 

7909 

945. 10 

0 

3 

0 

0 

0 

7910 

1226.35 

0 

7 

0 

1 

0 

7911 

1410.80 

1 

5 

0 

0 

0 

7912 

1286.15 

1 

5 

0 

0 

0 

8001 

1554.50 

3 

8 

0 

3 

0 

8002 

1902.45 

0 

31 

0 

2 

0 

8003 

2173.05 

0 

9 

0 

0 

0 

8004 

2633.55 

1 

2 

1 

1 

0 

8003 

2880. 15 

8 

9 

0 

0 

0 

8006 

3019.75 

1 

5 

1 

1 

0 

8007 

3423.05 

4 

12 

0 

1 

0 

8008 

3523.00 

0 

6 

0 

4 

0 

8009 

3193.60 

2 

7 

0 

2 

0 

8010 

4606.95 

2 

10 

1 

2 

2 

8011 

4150.80 

2 

9 

2 

0 

2 

8012 

3922.30 

1 

3 

0 

2 

1 

8101 

4288.90 

1 

4 

2 

0 

1 

8102 

4861.25 

4 

10 

0 

0 

3 

8103 

6093.80 

0 

19 

0 

0 

0 

8104 

6394.95 

3 

15 

1 

2 

1 

8103 

6384.50 

3 

14 

0 

4 

2 

8106 

7131.70 

1 

16 

1 

4 

0 

8107 

8212.45 

3 

18 

3 

1 

0 

8108 

1926.08 

12 

7 

0 

2 

1 

8109 

6686. 05 

7 

16 

4 

1 

0 

8110 

8748.00 

6 

14 

0 

2 

2 

8111 

8291.15 

3 

12 

1 

2 

2 

8112 

8884.20 

1 

16 

6 

1 

1 

8201 

9562.85 

7 

7 

3 

3 

5 

8202 

10150.75 

13 

26 

0 

1 

1 

8203 

10752.20 

15 

34 

3 

4 

4 

8204 

10941.95 

4 

35 

4 

1 

1 

8203 

11475.25 

9 

19 

4 

0 

1 

8206 

12829.65 

11 

39 

4 

3 

1 

8207 

12324.05 

4 

32 

3 

2 

1 

8208 

13628.70 

8 

57 

2 

4 

1 

8209 

12026.50 

9 

54 

0 

1 

1 

8210 

13626.40 

9 

36 

1 

9 

5 

8211 

14187.30 

3 

47 

4 

3 

1 

8212 

13188.40 

11 

64 

5 

4 

3 

8301 

13801.60 

13 

66 

5 

7 

2 

8302 

14361. 10 

18 

71 

2 

2 

14 

8303 

16156.15 

7 

97 

7 

4 

7 

8304 

15155.40 

10 

79 

8 

4 

6 

76 


Figure  B1 

Final  Periodogram  -  WUC:42GB0 
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Figure  B2 

Summary  o-f  Model  —  WUC:42GB0 
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Figure  B3 

Final  Periodogram  -  WUC: 14FB0 
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Figure  B4 

Summary  of  Model  -  WUC: 14FB0  El  Comp  Assy 
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Figure  B5 

Final  Periodogram  -  WUC:24EBA 
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Figure  B6 

Summary  o-f  Model  -  WUC:  24EBA  PTO  Sha-ft 
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Fi gure  B7 

Final  Periodogram  -  WUC:46AF0 
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Figure  B8 

Summary  o-f  Model  —  WUC:46AF0  Proport  Fuel  Flow 
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Figure  Cl 

Estimated  Inpulse  Response  -  WUC42GB0 

X  -  Prewhitened  Actual  Flying  Hours 
Y  -  Prewhitened  WUC:42GB0  Charger  AC  Battery 
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Summary  of  Model  -  WUC:42GB0  Charger  AC  Battery 
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Figure  C3 

Peri odog ram  -  WUC:42GB0  Charger  AC  Battery 

The  Estimated  Residuals 
Cumulative  Periodogram  .1  Probability  Limits 
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Figure  C4 
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Figure  C6 

Estimated  Impulse  Response  -  WUC:24EBA  PTO  Shaft 

X  -  Prewhitened  Actual  Flying  Hours 
Y  -  Prewhitened  WUC:24EBA  PTO  Shaft 


Figure  C7 

Sueeary  of  Model  -  WUC:24EBA  PTO  Shaft 
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Periodogram  -  WUC:24EBA  PTO  Shaft 

The  Estimated  Residuals 
Cumulative  Periodogram  .1  Probability  Limits 
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Figure  C9 


Histogram  -  WUC:24EBA  PTO  Shaft 
The  Estimated  Residuals 
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Figure  Cll 

Estimated  Impulse  Response  -  WUC: 46AF0  Proport  Fuel  Flow 

X  -  Prewhitened  Actual  Flying  Hours 
Y  —  Prewhitened  WUCS46AF0  Proport  Fuel  Flow 
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Figure  C12 

Suaaary  of  Modal  -  WUC:46AF0  Proport  Fuel  Flow 
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Figure  C15 

Cross  Cor r si st ions 
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Figure  C16 

Estimated  Impulse  Response  -  WUC: 46AN0  VI v  Shutoff  Mot  Opr 

X  -  Prewhitened  Actual  Flying  Hours 
Y  -  Prewhitened  WUC:46AN0  VI v  Shutoff  Mot  Opr 
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Figura  C17 

SuwMry  of  Modal  -  MUC:46AN0  VI v  Shutoff  Mot  Opr 
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Figure  C19 


Histogram  -  WUC: 46AN0 
The  Estimated 
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Figure  C20 


Cross  Correlations 
WUC: 46AN0  Proport  Fuel  Flow 
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